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Resumen 



Este trabajo esta dedicado al estudio de Teorias de Campo (QFT) en espacios 
difusos. 

Los espacios difusos son aproximaciones al algebra de funciones de un 
espacio continuo por medio de un algebra matricial finita. En el Kmite de 

matrices infinitamcntc grandes la aproximacion cs cxacta. 

Una caractcristica atractiva dc csta aproximacion es que mucstra dc una 
manera transparente como son preservadas las propiedades geometricas del 
espacio continuo. 

En el estudio del regimen no perturbativo de QFT los espacios difusos 
proveen una posible alternativa a la red como metodo de regularizacion. Esta 
tesis esta dividida en dos partes: 

1. Realizamos la simulacion Monte Carlo de la teoria en un espacio 
Euclideano de tres dimensiones. La regularizacion se compone por una 
esfera difusa de dos dimensiones, S"^, para las direcciones espaciales 
mas una red convencional para la direccion temporal. Se identifica el 
diagrama de phase de este modelo. Ademas de las fases desordenada 
y ordenada uniforme usuales encontramos una tercera fase de orde- 
namiento no uniforme. EUo indica la existencia del fenomeno conocido 
como mezclamiento UV-IR en el regimen de acoplamiento fuerte. 

2. Como segundo punto presentamos un analisis geometrico de una teoria 
escalar general en una esfera difusa de cuatro dimensiones, S%. Una 
aproximacion para es de un interes especial dado que S'^ es la subs- 
titucion natural de en estudios de QFT Euclideana. 
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Sin embargo la version difusa de S"^ no puede obtenerse mediante la cuan- 
tizacion del espacio clasico. 

El problema es rodeado definiendo una teoria escalar en un espacio mas 
grande, que es CP^ el cual es de dimension seis. Incluye grados de libertad 
de 5"^ mas otros que no lo son. Esos grados de libertad extras se eliminan 
dinamicamente mediante un metodo probabilistico. El analisis de las estruc- 
turas geometricas nos permite interpretar a este procedimiento como una 
reduccion de Kaluza-Klein de CP^ a S"^. 



Abstract 



This work is devoted to the study of Quantum Field Theories (QFT) on fuzzy 
spaces. 

Fuzzy spaces arc approximations to the algebra of functions of a continu- 
ous space by a finite matrix algebra. In the limit of infinitely large matrices 
the formulation is exact. 

An attractive feature of this approach is that it transparently shows how 
the geometrical properties of the continuous space are preserved. 

In the study of the non-perturbative regime of QFT, fuzzy spaces provide 
a possible alternative to the lattice as a regularisation method. The thesis is 
divided into two parts: 

1. We perform Monte Carlo simulations of a theory on a 3-dimensional 
Euclidean space. The regularisation consist of replacing space by a 
fuzzy 2-dimensional sphere, namely Sp, and Euclidean time by a con- 
ventional lattice. We identify the phase diagram of this model. In 
addition to the usual disordered and uniform ordered phases we find 
a third phase of non-uniform ordering. This indicates the existence of 
the phenomenon called UV-IR mixing in the strong coupling regime. 

2. Second we present a geometrical analysis of the scalar field theory on a 
4-dimensional fuzzy sphere, Sp. An approximation to S"^ is of special 
interest since S'^ is the natural replacement of in studies of Eu- 
clidean QFT. Nevertheless a fuzzy version of S'^ cannot be achieved by 
quantisation of the classical space. 



ix 



The problem is circumvented by defining a scalar theory on a larger space, 
CP^ which is 6-dimensional. It includes degrees of freedom related to S'^ plus 
others beyond S"^. Those extra degrees of freedom are dynamically suppressed 
through a probabihstic method. The analysis of the geometrical structures 
allows us to interpret this procedure as a Kaluza-Klein reduction of CP^ to 



Chapter 1 
Introduction 



Physics works best when there is a good interaction between experiment and 
theory. Unfortunately, for many of the interesting questions that arise, either 
it is impossible to perform the appropriate experiments or they are too costly. 
This is where the power of computer simulations plays an important role. In 
some sense, the computer simulation plays the role of the experiment. In 
modern particle physics, the experimental tools are large accelerators and 
the theories are typically quantum field theories. One of these is that of the 
strong interactions known as Quantum Chromodynamics. It is very difficult 
to extract some predictions from this theory as they fall in a non-perturbative 
regime and many physicists have resorted to computer simulations to extract 
the physical predictions. 

Similarly, many of the more speculative ideas emerging in physics involve 
strongly interacting field theories, some of these have novel features such 
as non-commutativity of the space-time coordinates. This type of structure 
is also suggested by string theory. The work of this thesis is dedicated to 
developing non-perturbative techniques adequate to these non-commutative 
theories and hopefully to string theory. 

Fuzzy spaces are included in the wider framework of non-commutative 
geometry. 

The idea of involving non-commutativity into Physics dated from the 
middle of the last century, nevertheless the substantial development has taken 
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place in the last few years. There are several reasons why studying non- 
commutative (NC) spaces has become so popular in the physics community. 
Although our interest in the study of fuzzy spaces is related to the study of 
Quantum Field Theories as we will discuss later, we would like to mention 
some other motivations for the study of NC geometry in Physics. 

Many interesting phenomena in Physics have been discovered by ex- 
tensions, therefore generalising commutative spaces into non-commutative 
spaces seems a natural extension. In this spirit, generalising commutative 
spaces to non-commutative spaces seems motivated. Non-commutativity 
can be incorporated into many branches of Physics like Gravitational The- 
ories, Condensed Matter Physics and Quantum Field Theories. The first 
attempts to involve NC theories in Quantum Gravity date from the last 
decade Fuzzy spaces can be found in String Theories (with D-Branes) 
under certain conditions — see Ref. [2]-|5. In Condensed Matter Physics it 
was found that the Quantum Hall Effect can be formulated in terms of non- 
commutative coordinates where a magnetic background field B is related to 
the non-commutative parameter [5]-[7j. 

One of the open problems in Field Theory is the existence of non well 
defined finite quantities: the divergences. The regularisation procedure mod- 
ifies the Field Theory to remove those divergences. The three well establish 
methods to the date are the dimensional regularisation, Pauli-Villars regu- 
larisation and the lattice regularisation [S]. The first two methods are for 
exclusive application at small coupling regimes. Regarding our motivation, 
we plan to test the feasibility of fuzzy spaces as a regularisation scheme in 
Quantum Field Theories. It should work, as the lattice procedure, at any 
regime. 

Suppose we want to study QFT through the path integral formalism on a 
given space. If we want to access the non-perturbative regime it is necessary 
to discretise the space in order to get a finite number of degrees of freedom. 
The standard method is to approximate the space by discrete points — a 
lattice — representing the space and then calculate the observables over that 
set of points. This simple idea has generated some of the most successful 
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theories in Physics, Lattice Field Theories, a review of them can be found in 
Ref. 0. 

On the lattice the continuous translational and rotational symmetry is 
explicitly broken. This certainly is a disadvantage in models where these 
symmetries are important. This takes us to search an alternative method 
that preserves these symmetries. 

Before proposing an alternative one has to ask which is the information 
necessary to describe an arbitrary space. The answer to this question was 
found by Connes and others in the context of non-commutative geometry [TU]. 
It is known that it is possible to reconstruct a manifolc^^ Ai if we have the 
algebra A of functions over Ai, a Hilbert space, Ti, and a differential operator 
able to specify the geometry (in Ref. [10] that operator is the Dirac operator, 
V, although for scalar theories as those studied in the present thesis it was 
conjectured in [H] that the Laplace-Beltrami operator A is enough to specify 
the geometry). Then, instead of discretising directly Ai by means of points 
we can discretise the triplet {A, H, V) and here the fuzzy spaces enter: they 
are essentially discretisations at algebraic level. If we want to obtain a finite 
number of degrees of freedom — namely, the coefficients in the expansion of 
a function in the algebra basis — the algebra has to be finite dimensional, 
i.e. a matrix algebra of dimension A^, Mat^, and as a consequence we have 
a non-commutative algebra. The elements in the algebra act on a finite 
dimensional version of the Hilbert space, T-Cn, and an appropriate version of 
T> is needed. In the limit — > oo (called the commutative limit) we have 
to recover Ai. 

Summarising, the fuzzy discretisation^ consist in the replacement: 

{A,n,V) {MatM,nN,VN). 

The term "fuzzy" originates from the following observation: since in the 
fuzzy space the coordinates will be matrices and they do not commute, this 

^Spaces are included in this more general notion of manifold. For our purposes we work 

only with spaces. 

^We denote it as fuzzification. 
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will mean in the Quantum Mechanics spirit that the notion of points does 
not exist, i.e. the space turns fuzzy. 

The seminal work on fuzzy spaces is due to Madore; in his work [12] a 
fuzzy approximation of a two-dimensional sphere is constructed. Since then 
there exists a large compendium of fuzzy literature, e.g. Refs. [13]- [IE]- Most 
of fuzzy spaces have been constructed based on the following observation: If 
we quantise a classical phase space we obtain a finite dimensional Hilbert 
space. This implies that the candidates to be quantised are manifolds of 
finite volume which have a symplectic structure. Co-adjoint orbits of Lie 
groups fall into this class. A didactic example of them are the complex 
projective spaces, CP", they are 2?7,-dimensional spaces that can be defined 
as SU{n) orbits. A discussion of its fuzzy version is given in Ref. |18J. The 
family of CP'^ is especially interesting since CP|, = Sp] Sp and Sp can be 
obtained form CP|, (see Ref. [T5]). 

Once we count with a fuzzy version of a space, the next step in our 
program is to define a Field Theory on it, e.g. Refs. [19] -[25]. Then we need 
to construct fuzzy versions of Laplacians, Dirac operators, etc. The solid 
mathematical background of the fuzzy approach makes it easy to identify 
such fuzzy versions. Of course one has to check whether the proposed theory 
reproduces the continuum theory. 

Field Theories on fuzzy spaces share a generic property of Field Theories 
on general non-commutative spaces called the UV-IR mixing [26] . The UV- 
IR mixing was originally discovered in perturbative calculations Refs. p7]- 
|29] . In non-commutative spaces we have two kinds of diagrams, those that 
reduce to commutative diagrams and those diagrams without a commutative 
counterpart, these are divergent at low momenta. 

Simulations on fuzzy spaces are a relatively recent topic, see Refs. [30] - 

m- 

In the present thesis we concentrate on the study of Scalar Field Theories 
on fuzzy spaces. We cover two important remarks in the fuzzy program: 

1. Test the feasibility as a discretisation method through a numerical sim- 
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ulation. 

2. Show that the fuzzy approach allows for a transparent geometrical anal- 
ysis of Field Theories. 

As a pilot study we describe a numerical study of the model on the 3 
dimensional Euclidean space S*^ ® T by means of the Metropolis algorithm. 
Our regularisation consists of 

• the fuzzy sphere Sp for the spatial coordinates 

• a conventional lattice with periodic boundary conditions for the time 
direction. 

The organisation of this first part of this thesis is the following: In chapter 
|4] we present a detailed description of both discretisation schemes, empha- 
sising the advantages of each method. In chapter [5] we present the charac- 
terisation of the phases in the model; we dedicate part of this chapter to 
the description of some technical aspects related to thermalisation problems 
in Monte Carlo simulations. In chapter [6] we identify the phase diagram of 
the model analysing the scaling behaviour of the critical lines. In chapter 
[7] we present the discussion of our results. The key point of this analysis is 
the behaviour of the triple point under different limits. It reveals that in 
the thermodynamic limit — > oo is it not possible to recover the Ising 
universality class due to the dominance of a phase that breaks the rotational 
symmetry spontaneously. We find that the UV-IR mixing predicted in the 
perturbative regime of the model appears in the strong coupling regime as 
well. There are ways to remove those divergences. In the context of the 
scalar field theory on the fuzzy sphere this is done by a suitable choice 
of the action [22] • The UV-IR mixing exists in the non-perturbative regime 
e.g. Refs. [35]-[38] where it was detected as a matrix or striped phase which 
has no counterpart in the commutative theory. 

As a second point in this thesis we present an analytical part. We study 
a Scalar Field Theory on Sp. 5*^ is a special 4-dimensional curved space, 
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taking its radius to infinity, we arrive at M^. is not a phase space, hence 
its "construction" involves some comphcations that are explained in chapter 
^ To solve this problem we allow S*^ to fluctuate into a larger space, this is 
CP^. We present a review of the construction of CP^ as SU{4) = Spin{6) 
orbit. Nevertheless for moving on 5*^, it is enough to preserve rotations in 5 
dimensions. We find that CP^ can also be constructed as a Spin{5) orbit, 
but demanding this less restrictive symmetry we construct a squashed CP^. 
For our purposes we demonstrate that this construction has more advantages 
since it allows us to identify CP^ as a fibre bundle over S*^ with S"^ as the fibre. 



We start chapter 10 with a short review of [39j. There, a Scalar Field Theory 
on an squashed CP^^ is defined. Then a penalisation method for all the non 
modes is introduced. This probabilistic method introduces an apparently 
"artificial" parameter h, such that h positive and large makes the non 
configurations improbable. Now we give an interpretation to this parameter 
through a geometrical analysis of the proposed model. Using coherent state 
techniques it is possible to extract the geometry. At the end we are able to 
"visualise" the penalisation method as a Kaluza Klein reduction of CP'^ to 
5*^. h is interpreted in terms of the radius of the fibre S"^. 



Part I 

Simulations of the Xcj)^ Model 
on the Space S% x S^i 
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Chapter 2 

Generalities of the method 



2.1 Path integrals and functional integrals 

The functional integral provides a powerful tool to study Quantum Field 
Theories. It can be thought of as a generalisation of the path integral for- 
malism in Quantum Mechanics introduced by R. P. Feynman in the late 40's 
(see e.g. [8] and |1Q]). The crucial idea behind the path integral is the su- 
perposition law. If we want to calculate the transition amplitude for going 
from an initial state at r' to a final one at time r", one has to consider a 
superposition of all possible paths. 

To state this in a mathematical form, let us suppose the initial state at r' 
to be denoted by {ipix')) and at time r" we have \ip{x")). Then the transition 
amplitude is given by 

{^l,{x")mr",r')mx')) (2.1) 

where 

U{r", t') = e-*^(-"-')/^ (2.2) 

Ti. is the Hamiltonian of the system, which we assume to be time independent. 
We start slicing the time interval [r', r"] into subintervals of duration 



Tj+j — Tj as in figure 2.1 



9 



10 



2.1. Path integrals and functional integrals 





To = T 

Figure 2.1: Schematic representation of the path. 



The time evolution operator can be broken into intervals: 



'l'H{T"-T')/h 



g- j;W(rjv-r]v-i+-rjv_i T2+T2--ti+ti-to) 



N 



(2.3) 



where r' = tq and r" = tn- 

2 

H has the general form H = Hq + V, where Hq = For e 



-in{T"~T')/h 



'iHoe/h -iVe/h' 



N 



(2.4) 



where we have used Trotter's formula. Eq. (2.4) holds if Tio and V are 



semibounded. The next trick is to insert between each term e ^'^o^/^e ^^'^/^ 
the set of complete states 



dxi\i){xi)){i){xi)\ = 1, 



(2.5) 



then eq. (2.1) can be written as the product of terms 



(V'(x")|Z^(T",r')|^(x')) 



dxN-i / dxN-2--- dx^{^{xM)\e-'^''"'e-''''"'\^P{xN-i)) 



\'4!{xN-2)) 



x{'4l{xN-l)\ 



x(V(xi)|e-'««^/V^^/''|V(xo)). 
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Suppose V depends on the position X and Ho depends on the momentum 
P. It is possible to demonstrate that (see e.g. Refs. 



(2.6) 



In eq. (2.6) we note that the argument in the exponential can be written as 



where we recognise a discrete version of Lagrangian in the interval [rj, Tj+i] 
times the duration of the interval e, i.e., the action in such an interval. Taking 



the product of the A^-terms of the type in eq. (2.6) and taking — > oo we 
arrive at 

{^{x')\U{T",r')mx'))= / [Z}x(r)] e^^W/^ (2.7) 

Jx{t')=x' 

where S [x] is the action 

S[x] = L{x,x)dT, (2.8) 

and L is the Lagrangian of the system. x{t) is a path that interpolates 
between x' and x", [Dx{t)] is the functional measure, therefore "J^ [Dx{t)]" 
denotes the integral over all paths between x{t') = x' and x(r") = x". 



An important remark is that from eq.(2.7) we can beautifully recover the 
least action principle noting that in the limit, h 0, the path of minimum 
action dominates the integral since the phase e*'^/^ of any path away from 
this fluctuates rapidly and different contributions cancel. 

The generalisation to quantum fields is a straightforward generalisation of 



eq. (2.7). But before introducing its expression we would like to remark that 
using path-integral methods it is common to give the action an imaginary 
time in order to simplify the calculations — the weight in the path integral 
is an exponential with real argument, which is easier to handle numerically 
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2.1. Path integrals and functional integrals 



— and then return to a real action at the end. This can be done if the 
Osterwalder-Schrader axioms hold (see Refs. 

Besides simplification purposes, in imaginary time the paths away from 
the classical path are exponentially suppressed. This makes the path integral 
to converge much better than the phase rotation. This is crucial for numerical 
studies since it allows us to have reliable results with a relatively modest 
statistics. 

There is a deeper consequence of considering an imaginary time: it allows 
us to establish a connection to Statistical Mechanics. 

2.1.1 From Euclidean time to real time 

In complex analysis, a branch of mathematics, analytic continuation is a 
technique to be used in the domain of definition of a given analytic function. 
We can apply such techniques here to go from real time r to the imaginary 
time t called the Euclidean time. (For a formal treatment see Ref. |H]). 

Imaginary time and spatial coordinates play equivalent roles. For exam- 
ple, in real time the D'Alembertian operator is given by: 

dr"^ dxl dxl dxl 
Under the Euclidean prescription we have 

dxQ dx\ dxl ^^V 

where we set xq = t. 

In Quantum Mechanics we consider the possible particle positions at each 
time, given by functions x{t) (or x{t) m d = 3), and the path integral inte- 
grates over all these functions, i.e. over all possible particle paths (with the 
given end-point). This reproduces the canonical Quantum Mechanics, but 
space and time are not treated in the same way. In field theory, one does treat 
them in the same manner and introduces a functions of any space-time point 
X = {x,t), which are denoted as fields. The simplest case is a neutral scalar 
field, where this field values are real, (plx) G M. The assignement of a field 
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value in each space-time point x is called a configuration, and it takes over 
the role of paths in Quantum Mechanics. Consequently, the functional inte- 
gral now runs over all field configurations, / [D(j)] . The Lagrangian L is now 
the integral of a Lagrangian density at each point x, (9^0(a;)), and the 

action is obtained by an integral over the space-time volume, S = J d'^xC. 
In imaginary time, the analog of eq. (2.7) for quantum fields is given bj|^ 



(0"|r(t",t')l0')= / [I^0(x,t)]e-^[^1. (2.9) 



In the Euclidean formulation the time evolution becomes a transfer ma- 
trix, T{t",t'), and in the case that H is time independent eq. ( |2.2 ) becomes 



g-wt_ Furthermore the quantum partition function Tr (^e~^^) becomes the 
functional integral over paths that are periodic in Euclidean time of period 
P, where (3 = with Boltzmann's constant. There is also a second 
interpretation of the resulting functional integral as a functional integral in 
statistical field theory. Here one considers the Euclidean action as the energy 
functional of an analog statistical mechanical system with ksT = 1. As is 
conventional in lattice field theory it is the latter analogue that will be used 



in this thesis. Then eq. (2.9) describes a statistical system in equilibrium. 



2.1.2 Expectation values 

The expectation values of an observable F, denoted (F), can be calculated 
as follows: 

{F) = ^l mF{<P)e-'^^\ (2.10) 

where 

Z = J [D0]e-^[^] (2.11) 

is the partition function. 



The integration in eq. (2.10) involves all the possibles configurations in 
the functional space. 



The problem is how to measure (or estimate) the value in (2.10). Here is 
where the importance sampling methods enters. The most popular approach 
^We set h=l. 
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2.1. Path integrals and functional integrals 



is the Monte Carlo Method (see appendix |A] for a brief description). The main 
idea of this method is that we can estimate (2.10 ) considering a representative 
sample of (independent) configurations. 

The way to produce the representative samples is through random moves 
to explore the search space. 

In this thesis we use a variant of the Monte Carlo method called the 
Metropolis algorithm [15] to estimate the expectation values of the observ- 
ables defined in chapter |4j 



Chapter 3 

A review of the 2 dimensional 
model on a fuzzy sphere 

We devote this section to a review of some aspects of the 2- dimensional 
model on a fuzzy sphere. We will discuss generic properties of fuzzy spaces 
by means of the most studied example: the fuzzy sphere. We will show that 
the fuzzy sphere can retain the exact rotational symmetry. 



3.1 The 2 dimensional Xcf)^ model 



A quite general scalar field theory on a 2 dimensional sphere is given by the 
action 



(0) := 



S2 



(3.1) 



where dVt = sinOdOdip, is a neutral scalar field on the sphere. It depends 
on the coordinates Xi{6,ip) which satisfy: 

3 



E 



(3.2) 



where R is the radius of the sphere. = ^^Cf, and Ci are the angular 
momentum operators. V [(f){x)] is the potential of the model. 



Eq. (3.2) describes 5*^ embedded in 



15 



16 



3.2. The fuzzy sphere 



3.2 The fuzzy sphere 

To obtain a fuzzy version of a continuous space we have to replace the algebra 
of the continuous space by a sequence of matrix algebras of dimension N, 
Mat AT. 

In the case of the fuzzy sphere, the permitted values of N are L + 1 where 
L is the largest angular momentum (the cutoff), which can take the values 
L = 0, 1, 2, ■ ■ ■ . The coordinates Xi are elements in the algebra of functions 
of S'^, C°°(5'^). They are replaced by the coordinate operators, X^, which are 
defined as = 2R ^^^_p where Li,i = 1,2,3, are the SU{2) generators in 
the N = {L + 1) -dimensional irreducible representation. 

The coordinate operators satisfy the constraint 

3 

J2xf = R'-l, (3.3) 

1=1 

which can be interpreted as a matrix equation for a sphere, the analog to 



eq. (3.2). Note that the operators X^ do not commute. 



2R 

[^.,^.] = ^e.,fc-^===X,. (3.4) 

Following the above prescription, the scalar field is represented by a her- 
mitian matrix $ of dimension A^. Just as in the standard case where can 
be expressed as a polynomial in the coordinates Xj, its fuzzy version $ can be 
written as a polynomial in the fuzzy coordinates. The differential operators 
Ci- are replaced by [Li, ■] and the integral over 5*^ is replaced by the trace. 
Summarising the above: 



Xi e C°°(5') 
Pix) e C°°(5') 
Ci(f)(x) 



S2 



{x)dn 



Xi e MatN, 
$ G MatN, 



AttR^ 



N 



Tr ($) 



C'. := J2 [L^, [L^, ■]] 



(3.5) 
(3.6) 
(3.7) 

(3.8) 
(3.9) 
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Note that G M implies that $ is hermitian and the choice of normah- 



sation in eq. (3.8) ensures that the integral of the unit function equals the 



trace of the unit matrix, i.e. 



N Js^ 



(3.10) 



Rotations on the fuzzy sphere are performed by the adjoint action of an 
element U of SU{2) in the dimension unitary irreducible representation. 
U has the general form U = e^^^\ The coordinate operators are then rotated 
as 

UXiU^ = RijXj, R G 50(3) (3.11) 



and the field transforms as 



(3.12) 



3.2.1 Limits of the fuzzy sphere 

Following |20], for the spatial part of our model (the fuzzy sphere) we have: 
• The commutative sphere limit S"^: 



N — y oo, R fixed. 



(3.13) 



The Moyal plane limit M.^ 
N — ^ 



oo, = — -— , 6 constant. 



The commutative flat limit 

N — > oo, R oc N-2^^-'K l>e>0. 



(3.14) 



(3.15) 



The limit given by eq. (3.13) arises naturally from the fact that N — > oo 
recovers C°°{S^). 



A short way to deduce eqs. (3.14)-(3.15) is the following 



18 



3.2. The fuzzy sphere 



Considering the north pole on the fuzzy sphere where X| ~ , we 



can re-scale the coordinate X| to -j^ 



relation (3.4) at the north pole as 

[^1,^2]-^ 



2R 



1 , and re-write the commutation 



2R^ X, 



, --X^i^i—-^. (3.16) 

We propose i? as a function in N. For the non-commutative plane we have: 



(3.17) 



then, comparing eq. (|3.16|) to (|3.17l) we obtain 



TV 



e. 



We define the exponent e in the relation 



R^ocN 



1 > e > 0. 



(3.18) 



If e = we have = const. For e > the commutator (3.16) vanishes if 
N — > oo. Note that in this limit we also require the commutator given by 



eq. (3.4) to vanish and this requirement is immediately satisfied for e > 0. 



3.2.2 The scalar action on the fuzzy sphere 

The next step is to define our field theory on the fuzzy sphere. Implement- 



ing the replacements given by eqs. (3.5)-(3.9) in eq. (3.1) we arrive at the 
following expression, 

47ri?2 '1 r2 



N 



1 n 



(3.19) 



Eq. (3.19) is valid for any potential ^ [$]. For testing purposes it is 
convenient to select a simple model. In Ref. ^32j the A</)^ model on a fuzzy 
sphere was studied, where the action is written as 



= Tr [a^C^(!? + 6<l>2 + c^^] 



(3.20) 



$ is a Hermitian matrix of size N. After a suitable rescaling, the parameters 
b and c become the mass squared and the self-coupling, respectively. 



In Ref. [32] was rescaled to fix a = 1. The model in eq. (3.20) was also 



studied in Ref. [30] but in terms of a different convention of parameters 

Air 



arR^ 



c = aXR . 



(3.21) 
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3.3 Numerical results on the two dimensional 
model 



The model in eq. (3.20) has been studied numerically by several authors 

We follow those results in Refs. ||32i, which are 



—see Refs. |30],|32]- 
summarised in figure 3.l|p 



2.5 



1.5 



0.5 



Disorder phase » 



N=2 I — I — I 
N=3 ^-x -j 
N=4 ; ^ ; ^ 
N=6 ^ B i 
N=8 
N=10 ; • 
cN"^=(bN"^'^)^/4 



Non-Uniform Order pfiase 



Uniform Order phase 



^1 

Triple point (0.80 ± 0.08,0.15 ± 0.05) 



10 



12 



14 



-b/N 



3/2 



Figure 3.1: Phase diagram obtained from Monte Carlo simulations of the 



model (3.20) in Ref. |32j. 



Figure 3.1 shows the existence of three phases: 



Disordered 



Non-uniform ordered 



• Uniform ordered 

^We thank the authors of Refs. [32j for their permission to reproduce the graph here. 
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3.3. Numerical results on the two dimensional model 



It is remarkable that the three coexistence lines collapse under the same 
re-scaling in for the axes 



vs. Two coexistence lines are: 



Disordered - Ordered uniform: 



-0.35ViVc. 



(3.22) 



Disordered - Non-uniform ordered: 



c 

iV2 



(6iV-3/2)i 



(3.23) 



The triple point is given by the intersection of the three critical lines 

(6^, c^) = (-0.15A^^/^ O.SA^^). (3.24) 



The transition in eq. (3.23) is only valid for large values of c. Therefore 



eq. (3.24) is not the intersection of eqs. (3.22)-(3.23). 



Although we discretised the spatial part of our three dimensional model 
by means of a fuzzy sphere, our simulations will show that the properties 



differ from those of the two dimensional model in eq. (3.20). In certain limits 
the X(f)^ theory on the fuzzy sphere can emerge as a limit of the 3-dimensional 



model in eq. (4.7) 



The model in eq. (3.20) depends on the parameters N, R, and A — or 
equivalently it depends on N, a, h and c, see eq .(3.21) — but it effectively 



depends only on two out of three parameters a, h and c. For the three 
dimensional model in addition to the parameters in the 2-dimensional model 
we have as parameters the number of lattice sites, Nt and the lattice spacing 
At. We will fix in section 4.4.3 Nt = N and the model will effectively depend 
on four parameters. 

The first question is if the three dimensional model has the phase of non- 
uniform ordering. We will see in chapter [5] that the answer to this question 
is yes. 



Chapter 4 



Description of the model 



In this Chapter we present the discretisation of the 3- dimensional model 
composed by a 2-dimensional sphere plus a Euclidean time direction. 

We first recall the results of chapter [3] for the discretisation of the continu- 
ous model. We will apply them to the 3-dimensional model. After performing 
the discretisation in the time direction we present the model to be studied 



by Monte Carlo techniques. In section |4.4| we present the observables and a 
brief description of their meaning. 



4.1 Regular isat ion of the action 



It is convenient to consider the Euclidean version of the model. 

As it was remark in chapter |2} the main advantage of working in this 
formalism is that it allows to establish a connection to Statistical Physics 
and the functional integral converges with a relatively modest statistics. 

The model to regularise is 

,2 



dt 

51 ^52 



£2 



{x,t) + "^4>\x,t) + ^cl)\x,t) 



(4.1) 

X, t) is a neutral scalar field on the sphere. It depends on time (euclidean) 



and on the coordinates Xi{6, (p) satisfying eq. (3.2), where R is the radius of 
the sphere. 
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4.1. Regularisation of the action 



We consider the integral over Euchdean time on a compact version, Si, 
which has circumference T. 

First we will explain how to discretise the spatial directions and then we 
will perform the discretisation in time direction. 



4.1.1 Discretising by a fuzzy sphere 



Let us consider the spatial part of the action given in eq. (4.1) 
'1 



52 



Tfl A 

x,t) + —4>''{x,t) + -<p\x,t) 



R^dn. 

(4.2) 



Implementing the replacements given by eqs. (3.5)-(3.9) in eq. (4.2) we 
arrive at 



^,-dnm + '^^\t) + \^^(t)] . (4.3) 



Then, the action (4.1) "discretised" in the spatial directions is: 



2^^ W - <^> it) + ^'^\t) + -1>\t) . (4.4) 



The model given by eq. (4.4) has the exact rotation symmetry of model (4.1 ) 



since any rotation on the sphere is allowed and the action (4.4) is invariant 



under uniform rotations given by eq. (3.12) 



4.1.2 Discretisation of the time direction 



To discretise the time direction we take a set of Nt equidistant points, then 
T = NtAt. 



The changes to implement in eq. (4.1) are: 

dt - 



Nt 



'Si 

dt(j){x,t) 



t=i 



(x,t + At) - (j){x,t) 
At ' 



(4.5) 
(4.6) 
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We arrive at 



t=i ^ 



(4.7) 



One configuration $ corresponds to a set of matrices {$(t)}, for t = 1, . . . , Nt. 



Alternatively we can write down eq. (4.7) in terms of the constants A, D, B and 



C defined in eqs. (4.8|-(4.11 ): 



A 
D 
B 
C 



2TTAt 

IT' 

2ttR^ 
NAt' 
2TTR^m'^At 
N ' 
TTR^XAt 



N 



(4.8) 
(4.9) 
(4.10) 
(4.11) 



Then the action reads: 

Aft 



= ^TY^A$(t)£2$(t) +D($(t + At) - $(t))VB$2(t) + c$^(t) . (4.12) 



t=i 



4.2 Decomposition of the field 

As we mentioned in the previous section, we are representing a configuration of 
the field in our model by a set of matrices {^(t)}, for t = 1, . . . , Nt. Every element 
in this set can be expanded in the polarisation tensor basis 

N-l I 

m = E E ^imit)yim, (4.13) 

1=0 m=-l 

where Qm(t) are A^^ coefficients. The polarisation tensors Yi^ are N x N matri- 
ces that are the analog of the spherical harmonics, Yijn{9,{p). Details about the 
polarisation tensors are presented in appendix |B| 

At the end, the quantities of interest can be expressed as expectation values 
or averages over the configurations. The expectation value of the observable -F($) 



was defined in eq. (2.101: 



-sm 
z 



(4.14) 
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4.2. Decomposition of the field 



where 2 = J [D^] e ^^^^ is the partition function. Quantities of interest wih be 



mmt')) 



(4.15) 
(4.16) 



where eq. (4.15) is a condensate and eq. (4.16) is a correlation function. These can 



be mapped to standard correlation functions by replacing the Yim by Yim{9,ip)- 
We can reduce the expressions ( |4.15 )-(4.16 ) to combinations of the expectation 



values of the coefficients cim{t) introduced in eq. (4.13) 



{Clmit)), 
{cim{t)cVm'{t')). 



(4.17) 
(4.18) 



Now it is convenient to compute the quantities (4.17)-(4.18) after a Fourier trans- 
form in (Euclidean) time. 

Following Ref . |29] , the complete Fourier decomposition of the field is given by 



Nt-l 
l,m k=0 



(4.19) 



where 



^ t ^ 



(4.20) 



Clm{t) 



It will sometimes prove convenient to define 

1 



27rNt 



e $(t). 



(4.21) 



In this space the correlator (4.18) is diagonal, 



(4.22) 



Gim{k) is the Green function in momentum space. Here the term "momentum 
space" is used to include both angular momentum {l,m) and frequency k. 
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4.3 The different limits 



For all our simulations we are interested in taking the thermodynamic limit N — > 
oo. 

For the time direction, we are interested on taking Nt — > oo. 
We now set to At = 1. Now, the tricky part is how to relate the parameter R, 
the radius of the spheres, and A^, the dimension of the matrices. 



For the limits of the spatial part of our model we follow section 3.2.1 



4.4 Definitions of the observables 



We define the field averaged over the time lattice as: 

The average over the time lattice of the coefficients ciui are: 



Clr 



(4.23) 



(4.24) 



This picks out the zero frequency component of ^, i.e. eqs. (4.23)-(4.24) are 



particular cases of the equations (4.19 1-(4.21 ) when k = 



Some particular cases in eq. (4.24) are 

\/47r. 



Coo 

Clm 



N 



Tr$, 



N 



(4.25) 
(4.26) 



where li,m are given in eqs. (B.12 )-(B.14 ) of appendix B 



4.4.1 Order parameters 

We want to measure the contributions of different modes to the configuration 
For this purpose we need a control parameter. This turns out to be the sum 

|c;mP) this quantity was called the full power of the field in Ref. [30] and it 
represents the norm of the field it can be calculate as : 

^lir-=E\^im\' = %TT{^'). (4.27) 

l,m 
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4.4. Definitions of tlie observables 



Although (95^;;) cannot play the role of an order parameter, we will show that 
it is useful to localise the region where the phases split into disordered and ordered. 



We expect {^p] 



all I 



in the disordered phase and (y?^;;) 



S> in the ordered phase. 



To distinguish the contributions from the different modes to eq. (4.27 ) we define 
the quantity: 



'fl 



We can re- write eq. (4.27) in terms of the quantities in eq. (4.281 



2 

Vail 



(4.28) 



(4.29) 



In the disordered phase we expect {(pi) ~ for all /. 

Studying the contributions of the different modes to (92^;;) can provide more 
information about the phases. If ((/?/) ^ for / > it indicates that the rotational 
symmetry is broken. In our simulations we measure quantities related to the lowest 
modes: the zero mode for I = and the first mode for Z = 1 as representative of 
those modes where the rotational symmetry is broken. 



Choosing the particular case / = in eq. (4.28) we have 



ipo := I Cool- (4.30) 

For m? < 0, if the contribution of the fuzzy kinetic term to the action is not 
negligible we can expect the kinetic term to select the zero mode as the leading 
one, {(Pq) = i^Paii)- As a consequence {ipo) ^ in the uniform ordered phase; (tpi) 
is expected to be close to zero in the disordered phase. 
Its corresponding susceptibility is defined as: 



Xo := ivl) - ivof. (4.31) 

As the contribution of the kinetic term to the action reduces compared to the 
potential contribution we can expect the system can undergo the condensation of 
higher modes. Let us consider the p-wave contribution to ^, i.e. the contribution 
of the I = 1 mode. Using cim, m = 1,0, —1 we introduce a 3-dimensional vector, 

( ^ 
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With this vector we can define the order parameter (fi, as a particular case 



Z = 1 in eq. (4.28) we have 



and its susceptibility, xi- 



\ 



1 

E 

m=—l 



Cl, 



12 .- 



(4.32) 



XI := {ifi^) - {ifif. (4.33) 

Following Ref. [30] the ordered non-uniform phase is then characterised by 
{(fi) ^ 0. Note, however, that due to fluctuations we will always have {ipf) > 0, 
so we have to specify how large it has to be. We will give more details of how to 
characterise this phase in the next section. 



We can include contributions of the remaining modes generalising (4.32) and 



(4.33 ). In practice the study of the first two modes should be enough to understand 



(4.34) 



the behaviour of the system. 

4.4.2 Energy and specific heat 

The internal energy is defined as: 

E{m^X) := (5), 

and the specific heat takes the form 

C{m^,X) := {S^) - {Sf. (4.35) 

These terms correspond to the usual definitions E{m?, A) = and C{m^, A) = 

where Z is the partition functionj^ 

We separate the action ( 4.12[ ) into its four contributions: 

t 

2 



dE 

80 



Si m 

52[$] = D^TY(^$(t + l) -$(t))', 

t 

= B^Tr(^'^\t)), 

t 

S,[^ = cY,T^rU\t)), 



(4.36) 
(4.37) 
(4.38) 
(4.39) 



^[3 is proportional to the inverse of the temperature T, i.e. /3 — fr^, where fcs is the 
Boltzmann constant. 
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4.4. Definitions of tlie observables 



where A, B, C,D were defined in eqs. (4.8)-(4.11 ). 



Tlie corresponding expectation values to eqs. (4.36 )-(4.39 ) are 



E2{m^A) 



(Si) 
{S2) 
{S3) 
{S4) 



(4.40) 
(4.41) 
(4.42) 
(4.43) 



4.4.3 Dimensionless parameters. 



Eq. (4.7) is written in terms of tlie following parameters: a general temporal 
lattice spacing At , the radius of the sphere R, the dimension of the matrices 
A^, the number of points in the lattice A^t,the mass squared and the self- 
coupling A. 

In order to simplify the simulations, we use the freedom to re-scale the field 



<I> to fix the value of one of the constants given in eqs. (4.8)-(4.11 ). For our 



simulations we fix ^ = 27r — see chapter [7] for more details. We defined the 
dimensionless parameters: 



m 



R 

At' 
AtA. 



(4.44) 
(4.45) 
(4.46) 



In all our simulations for the 3-dimensional model Nt was taken equal to N. 



Chapter 5 



Description of the different 
phases in the model 

In this chapter we characterise the different phases present in this model. Accord- 
ing to the values of A relative to a critical value Xx we will see that we can divide 
the space of parameters into two regions. In both cases we can subdivide according 
to values of m?: 

1. At > A > 0. 

(a) For < we have a uniform ordering (Ising type). 

(b) For > ml we have the disordered phase. 

2. A > At. 

(a) For < we have a non-uniform ordering. 

(b) For > we have the disordered phase. 

5.1 Behaviour of the system for Ay > A > 

In the previous chapter we defined the observable if'^n called the full power of the 
field. It provides a control parameter since it represents a norm of the matrix 
We will sec that {^'^n) ~ defines a disordered phase while {(f'^n) ^ defines 
a kind of ordered regime. But to describe the type of ordering it is necessary to 
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5.1. Behaviour of the system for At > A > 



study the contribution from separate modes to (v^^iz)- figure 5.1 we present 



I — \ — I 



a typical case for > A where we show the partial contributions from the zero 
and first mode to {^p'^n)- 
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Figure 5.1: {ipln), {iff) and {^f) vs. at A = 0.17, i? = 4, iV = 12. 



From figure 5.1 we can observe that for m? > —0.1: ip'^ 



all 



0, m2 < -0.1: 



^aii ^ fo > ^1 SO the dominant mode turns out to be the zero mode. 
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Figure 5.2: {(po) and (v?i) vs. m2 at A = 0.17, ^ = 4, = 12. 



Figure 5.2 shows the order parameters {(fo) and {(fi). In order to determine 
precisely where the phase transition occurs the standard way is to search the 
maximum in the susceptibility, in this case xo since the zero mode is relevant for 
this phase transition. 



5 Description of the different phases in the model 



31 



-2 







0.018 
0.016 
0.014 
0.012 
0.01 
0.008 
0.006 
0.004 
0.002 




-O.i 



-0.6 











1 

Xo 1 — 1 — 1 









\ 


\ 




























































+ + 


















■H 











-0.4 



0.2 

"2 



0.2 



0.4 



Figure 5.3: The susceptibilities xo and xi, in eqs. (4.31) and (4.33), at A 
0.17, i? = 4, = 12 



Xo peaks at rn? = —0.1 and in figure 5.3 we can observe that the susceptibihty 



associated to the first mode reveals a small response too. 



Figure 5.4 shows the internal energy for the same parameters as in figures 
[53D 
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Figure 5.4: Internal energy density see eq. (4.34), and its partial contri- 
butions, given in eqs. (|4.40[)-(|4.43[), at A = 0.17, ^ = 4, iV = 12. 
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5.1. Behaviour of the system for At > A > 



We observe in figure 5.4 that the leading contribution for m? > —0.1 to the 



internal energy E, is the one that comes from the kinetic fuzzy term in eq. (4.40). 
Note that the phase transition occurs at m? where the potential contributions £"3 
and -E'4 deviate from zero. 



Figure 5.5 shows an archetypical behaviour of the specific heat for > A. 
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Figure 5.5: Specific heat per volume at A = 0.17, i? = 4, = 12. It 



follows from eq. (3.10) that for this values of the parameters the volume is 



the constant AttR'^ x Nt = 19277. 

The specific heat provides an alternative criterion to the susceptibilities to deter- 
mine where the phase transition occurs. It provides information about the order 
of the phase transition. We prefer at this point to follow the specific heat criterion 
because it is a more universal quantity that does not distinguish the dominant 
mode. Since we expect that as we increase A the dominant modes are higher than 
the zero and first mode, we cannot ensure that in that region the susceptibilities 
related to such modes xo and xi give a reliable prediction to the critical point. 
If we follow the criteria of the susceptibilities we have to take into account which 
mode is the dominant one. Both criteria are theoretically supposed to detect the 
same phase transition at the same parameters. For A < At where the dominant 
mode in the ordered phase is the zero mode, this is confirmed comparing figure 



5.3 to figure 5.5 since the susceptibility of the dominant mode, xo: and the spe- 
cific heat peaks around the same value in m'^. We can also observe in figure 5.5 
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that there is a smaller response in the susceptibility of the first mode, xi- For 
A < At we will see that this situation is different since the susceptibilities of the 
non-dominant modes do not peak at the phase transition. 

We call the value of where the specific heat peaks m^, and for figure 5.5 
= — 0.12ib0.02j^ In the case that for > the kinetic term is not leading, it 
appears as a small shift between both peaks. This happens when R is big enough 

the error 



to have ^ small, At > A. Another observation is that in figure 



5.3 



bars are smaller than in the case of the specific heat from figure 5.5 The reason 
is that in general more statistics is necessary for the specific heat than for the 
susceptibilities. This phase transition is of second order as it is shown in figure 
[53] 



5.1.1 Thermalisation with respect to the observables 

In this section we want to present the typical behaviour of the observables for 
> and < m^. 



First we discuss some aspects of the thermalisation procedure. We define the 
thermalisation time as the number of Monte Carlo steps necessary for an observable 
to stabilise around one value independently of the starting conditions. If > 
we are in the disordered phase that is characterised by the property that the 



coefficients c;^ in eq. (4.24) are in average near to zero. We assume that the 
thermalisation of the coefficients c/^ is similar, and we check if the coefficient cqo 
thermalises. Figure 5.6 shows the thermalisation of the action and the coefficienlj^ 



Cqo for the point = in the figures |5.lf |5.5[ 



^For practical purposes we had to estimate the error by referring to the spacing of the 
values that we simulated. 

^During the run the values stored were Cqo and from them we can trivially calculate 
ifQ. We prefer to present the histories and histogram of Cqo rather than ipo in order to 
check if the samples are symmetric under Cqo — > ^cqq- 
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5.1. Behaviour of the system for At > A > 
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Figure 5.6: Example of thermahsation of the action and the coefficient Cqo 
at m2 = 0, A = 0.17, i? = 4, iV = 12. 



We chose a hot start |^ for figure 5.6 . We simulate the same parameters with 
a cold start 1^ and we obtained results in agreement within the statistical errors. 



The thermahsation time for cqo in figure 5.6 is estimated to 1500 Monte Carlo 
steps. 



After the thermahsation we begin the measurement procedure. Figure 5.7 
shows the histogram^ for the observables and parameters used in the simulations 
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Figure 5.7: Histogram of the of the action and the coefficient Cqo at = 0, 



A = 0.17, R 



4, N 



12. 



If > we are in the ordered regime, for A < At we have uniform ordering 



characterised by the property that in the expansion (4.23) all coefficients cim for 



^The starting configuration is a vector of hermitian matrices filled in with random 

numbers. For more details see appendix [X| 

^The starting configuration is a vector of matrices proportional to the unit. 
^The area is normalised to 1. The number of bins is 500. 
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I > average zero. Again it is sufficient to clieck if tlie coefficient cqo tliermalise. 



Figure 5.8 sliows tlie tliermalisation of tlie action and tlie parameter Cqo for tlie 



point 
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Figure 5.8: Example of the thermalisation of the action and the coefficient 
Coo at m2 = -0.2, A = 0.17, = 4, iV = 12. 



We chose a hot start for figure 5.8 After 1500 Monte Carlo steps the value 
of Cqo oscillates around 0.9 and the energy per unit of volume fluctuates around 
0.43. We also simulated at the same parameters with a cold start . We obtained 
the same results within the statistical error but the thermalisation time decreases 



by more than 50% as we can observe in figure 5.9 
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Figure 5.9: Example of thermalisation with cold starting conditions of the 
action and the coefficient Coo at = —0.2, A = 0.17, .R = 4, = 12. 



Finally, after the thermalisation procedure we measure the expectation value of 
the observables. We present the histograms for the same parameters as in figure 
[531 
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Figure 5.10: Histograms of the action and the coefficient Cqo at = —0.2, 
A = 0.17, R = 4:, N = 12. 



The peaks in th e prob ability distribution of cqo are approximately located at 
the peaks are approximately located at cqo 



Infi 



2 

47rx0.2 
12x0.17 



gure 



5.10 



±0.95 and 



1.12. As we move forward for a more negative this prediction is 



more accurate. This is shown in figure 5.11 
70 




coo/\/47r X 2/(12 x 0^17) 
Figure 5.11: Histogram of coefficient Cqo at = 

= 12. We re-scale the x-axis by the factor 

takes the value 3.54. 



-2, A = 0.17, R = 4, 



4-7r|m2| 
NX 



which in this case 



5.2 Behaviour of the system for A > At 

For this region of the phase diagram we observe two phases: for > we 
have the disordered phase characterise by {^'^n) ~ 0, for < we have the 
ordered phase characterised by {f'^n) > 0. For <C there are thermalisation 



problems, we will discuss these difficulties in section 5.2.1 
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In Figure 5.12 we present the partial contributions from the zero and first 



mode in (iflu) 
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Figure 5.12: {ipl^), {<^l) and {y:>l) vs. m2 at A = 1.25, i? = 8, = 12. 



We observe in figure 5.12 {ip'^ 



all/ 



for > —0.3 while for < —0.3 we have 



i^aii) ^ *-* ^^"^ —0.6 > > —0.1 we observe {faii'^) 








Figure 5.13: {(fo) and {(fi) vs. m2 at A = 1.25, R = 8, N = 12. 



For those values of where < 
the quantities ((/9q) and {^pf) in figure 



0.6 we can observe a kind of irregularity in 
and in ipo and <pi in figure 



5.12 



5.13 



smce 



they do not grow monotonously. We will come back to this point at section 5.2.1 



For the moment we focus on the region where the observables behave smoothly, 



i.e. > —0.6. 
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5.2. Behaviour of the system for A > At 
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Figure 5.14: Susceptibihties xo and xi at A = 1.25, i? = 8, = 12. 



In figure 5.14 we can observe that xi indicates a phase transition for ~ 
—0.35 while xo cannot detect it since xo keeps growing as tit? decreases. We 
conjecture that xo should peak for some value of < —0.6. 
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Figure 5.15: Internal energy E of eq. (4.34) and its partial contributions 



eqs. (4.40)-(4.43) at A = 1.25, i? = 8, iV = 12. 



The specific heat in figure 5.16 indicates a phase transition at = —0.37 it 
0.02. 
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Figure 5.16: Specific heat at A = 1.25, R = 8, N = 12. 



As we can observe comparing figure 5.14 to figure 5.16, there is a small 
difference in the critical value of predicted by the susceptibility of the zero 
mode xo the one given by the specific heat. We conjecture that this is due 
to a finite volume effect, nevertheless both criteria are qualitatively the same. For 
more details see appendix [Cj 



5.2.1 Thermalisation problems 



In this section we want to sketch the thermalisation problems. 

The program was designed to perform an arbitrary number of independenlj^ 
simulations in every run, Usim- In the case of figure 5.17 we performed ten 
independent simulations, the first three of them with a cold start while for the last 
seven simulations we chose a hot start. 



^If we choose a hot start, they have different starting configurations. 
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Figure 5.17: History of the action and the coefficient cqo for different starting 
conditions at 



-0.66, A = 1.25, R = 8. N = 12. 



We can observe in figure 5.17 that the expectation value of the energy and the 
coefficient cqo depend on the starting conditions. The interpretation of this phe- 
nomenon is that the effective potential of the system has several local minima with 
barriers large enough to suppress tunnelling between them. For different starting 
conditions the system gets trapped in one of those minima. As a consequence we 
have thermalisation problems (or practical ergocidity problems in the algorithm). 
They should disappear at an infinitely large Monte Carlo time, Tmc- In figure 
5.17| we have a new simulation every 1, 500, 000 Monte Carlo steps from which the 
ffi'st 500, 000 steps where taken as thermalisation time. 

The large error bars in the observables at = —0.66 in figures 



5.12 



and 



5.13 can be explained because the different simulations give different results. For 
cold starts the action of the system oscillates around 0.05 - this is the value of the 
energy in the absolute minima - and the trace cqo of the sampled configurations 
fluctuate around 0.17; for hot starts the energy of the system oscillates around 0.18 
and the trace cqo of the sampled configurations fiuctuate around zero. 
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5.3 Estimating the maximal number of min- 
ima 



It is possible to observe, from the Monte Carlo time evolution of the observable 



Coo in ec[- (4.25) that there are no thermalisation problems and that its probability 



distribution has several peaks. 
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Figure 5.19: History of the coefficient Coo for 
N = 12. 



6000 



0.3 at A = 0.75 R = 16, 



The history corresponding of the action for the same parameters in figure 
5.19 is the following: 
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5.3. Estimating the maximal number of minima 
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Figure 5.20: History of the action for 
= 12 



-0.30 at A = 0.75 R = 16, 



In this case the fluctuations were large enough to jump from one minimum to 
another. Nevertheless the configurations sampled belong to different subspaces of 
the spaces of configurations characterised by the different values of cqo it is shown 



in figure 5.19 



The histogram corresponding to figure 5.19 is presented in figure 5.21 
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Figure 5.21: Histogram of the coefficient Cqo for = 
^ = 16, = 12. We divided the x-axis by the factor - 
to emphasise that the peaks are around integer values. 



0.30 at A = 0.75 
^ 0.0045 



The estimated value for the triple point for N = 12, i? = 16 is (Ar,?7i^) 



(0.25, —0.062), the parameters simulated for figures 5.19 5.20 are in the region 
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A > At (i.e. in the ordered non-uniform phase). For this region the leading 
contribution comes from configurations characterised by a angular momentum 



Z > as we can check from the expansion of the norm of the field eq. (4.27): 
V'aii ~ V'o + V^i + ■ ■ ■ parameters simulated for figure 

tions to {iflii) split as follows {ipli^) ~ 0.30 ~ 0.014 + 0.23 + 



5.19 



the contribu- 



The several maxima in figure 5.21 reflect the existence of several minima 
in the effective action. If A ^ Ay it is expected that the potential provides the 
leading contribution to the action, turning into a pure potential model. Then, if 
A is large enough, for this region of parameters the minima in the effective action 
are given by the minima in the potential. 



5.4 The equilibrium configurations 

The maximal number of local minima for our model can be obtained from the 
maximal number of local minima for the two-dimensional model. In Ref. |30) it 
was conjectured that the minima in the 2-dimensional pure potential model are 
given by the disjoint orbits: 

,2 



o„ 



{^U^ (In e l2s+i~nU\U e U{2s + l)/[U{n) X U{2s + 1 - n)])}, (5.1) 



A 



where n < s + ^ and In are n x n matrices (2s -|- 1 = A^), i.e. we have N disjoint 
orbits. 



In our case that we have Nt lattice points. Note that in the expression (5.1 ) the 

o2 



constant 



A 



is in terms of the parameters of the 2-dimensional — see eqs. (3.20) 



(3.21). For large values of A we can establish their "equivalence" in the three 



dimensional model via eq. (6.51) — see section 6.5.2 We conjecture that the min- 



ima in the 3-dimensional model are at 



m)} 



Nt 
t=l' 



iVA 



■U^AU 



(5.2) 



where A is a diagonal matrix with eigenvalues 1 or — 1. 



To prove the expression (5.2) we analyse cqo 



The coefficient cqo ~see eq. (4.25)- can be expressed as 



Coo 



NNt 



Tr 



Nt 



.t=i 



(5.3) 
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5.4. The equilibrium configurations 



From expression (5.2) we have NNt + 1 possible values 

-NtN, -NNt + 2,..., NNt -2, NNt, 




Nt 



(5.4) 



t=i 



then 



v/47r|m2|/(ArA) 



NN 



X Coo = Tr 



Nt 



.t=i 



takes integer values as in figure 5.21 where N x Nt — 2k 
The maximal number of minima is x A'^ + 1 . For N 



(5.5) 

-12,-6,0,6,12. 

12 = A'f as in figure 



5.21 the maximal number of peaks is 145 but we just observe 5 of them. 



We can discuss intuitively why we cannot observe the maximal number of 
peaks. There is a single way to obtain the value -v/ ^ X^t^i Tr($(t))A^ x Nt in 



eq. (5.4) -all the matrices $(t), t = 1, • • • ,Nt should be proportional to the identity. 
To get in eq. ( [5^ the value NNt-2 we should have ^>(to) oc U'>Diag{l, 1, • • • ,1, -1)U, 
$(t) oc WlU for t 7^ to- Since X^^x '^('^(^)) invariant under interchange of the 
lattice points and interchange of the eigenvalues on each matrix Tr(<I>(t)) there are 



NNt ways to obtain the value NNt — 2 in eq. (5.4). This give us an rough notion 



of the volume of the potential minima. Those when eq. (5.4) are around have a 
larger volume than those where Tr($(t)) = NNt, they are the most probable to 
fall in a simulation with hot starting conditions. 

Since the number of local minima is larger the measurement problems are more 
severe than in the 2-dimensional case, see e.g. Ref |30]-|32]. The tunnelling between 
different minima depends on the size of its potential barrier and the size of the 
fluctuations. 

In general we can diagonalise only one of the matrices {^(tj)}^^. This 
is represented in figure 5.22 where the matrix ^{tj) is diagonalised though a 
rotation U £ SU{2) to Aj. For {^{ti)}^^, i ^ j the matrices are not diagonal. 

^{tj) oc Aj it can be map to the continuum to f{0, if) 



N-l 

E 

1=0 



N-l 



<!>{tj) OC Aj = fiYio f{9,ip) = fiYioie,^). 



(5.6) 



1=0 



This is schematically shown in figure 5.22 with the sign^"ib" in red on ^{tj). 



'^Yio{0, V3) > for f > 61 > and ¥10(6, (/?) < for tt > > f . 
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Figure 5.22: Schematic view of the equihbrium configuration. 



Note that the model in eq. (4.7) can be brought to one of diagonal matrices if 
the contribution of fuzzy kinetics term in eq. (4.36) is negligible (see Ref. |48]-|49]). 



We have the picture it figure 5.23 




Figure 5.23: Schematic view of the equihbrium configuration when the fuzzy 
kinetic term is neghgible. 



Chapter 6 

The scahng behaviour 



In this chapter we focus on finding out the dependence of the transition curves 
and triple point on the parameters of the system. 



6.1 Phase transition disordered to ordered- 
uniform 

First we explore the transition curve for = 12, R = 8 and we consider several 
values for A. 



I - II : N =12, R-- 

-0.410A 
-0.426A1050 

I: Disordered 




Figure 6.1: Transition curve from the disordered to ordered-uniform phase 
for N = 12. 
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6.1. Phase transition disordered to ordered-uniform 



We observe in figure 6.1 that the transition curve is a line that crosses the origin, 
given by the equation 

m2 = -(0.410 ±0.007)A. (6.1) 
In addition we tried a fit of the form m? = aX^ and we found the expression 

m2 = -(0.426 ± 0.021)A(^-°5°^°-°6°). (6.2) 



Eqs. (6.1 1 and (6.2) are in agreement. We choose a linear fit for our transition 
curves. 

For general values of N and R we propose a transition line between the disor- 
dered phase and the ordered uniform phase of the form 

m^ = h{N,R)\. (6.3) 
To identify the slope fi{N, R) our strategy is the following: 

1. First we extract the dependence on varying N and keeping R fixed. What 
we would expect is of the form 

m2 = fi{R)N^^X. (6.4) 

2. Then, if 5i does not depend on R we extract the dependence on R proposing 

h{R) = const. R^\ (6.5) 



3. Note that eq.(6.3) would still be valid if both sides of the expression are 



multiplied by a common factor. We will use this factor to stabilise the triple 



point to a fixed value. This procedure will be explained in section 6.3 



The slope /i ( A^, R) could have a more complicated dependence on N or R, 
nevertheless we will prove the choice fi{N,R) = const. N^^R^'^ is a good ansatz. 

Step 1 

Now, to find the collapse on N , we keep R = 8 and we consider N = 8, 12, 16, 23, 33. 
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Figure 6.2: Transition curve from the ordered- uniform phase to the ordered 
non- uniform phase for R = 8. 



In figure 6.2 we re-scale the j;-axis by a factor of N^'^^. We will see in section 



6.3 this factor stabilises the value of the triple point in N. 



The equation of the fit of figure 6.2 is: 



-0.0817V°-^^A. (6.6) 
The exponent on N, 5i = 0.64 is optimal for this fit; to compare with another 



exponent we show the same data on the next figure 6.3 for di = 0.5: 
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Figure 6.3: Transition curve from the ordered- uniform phase to the ordered 
non- uniform phase for R = 8 for 6i = 0.5. 
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6.1. Phase transition disordered to ordered-uniform 



As we observed comparing figure 6.2 to figure 6.3 the fit for the exponent 
5i = 0.64 is better and we estimate the error on 6i to the value 0.2 for R = 8, this 
is 6i{R = 8) =0.64 ±0.2. 



The next case R = Am figure 6.4 seems to confirm the value for the exponent, 
5i{R = 4) = 0.64 ±0.2. 



— 1 — 1 

iV = 33; S = 4 

iV = 23;S = 4 >— A-j. 

TV = 16; B = 4 :- - • - ■: 

iV = 12;_R = 4 r--Q --. 
Ar = 9;fi = 4 
Ar = 8;i? = 4 ^ 
-0.127N°-^*X 




I: Disordered 



Figure 6.4: Transition curve between the disordered phase and the ordered- 
uniform phase for ^ = 4. 



We illustrate the case i? = 16 in figures 6.5 6.6 
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Figure 6.5: Transition curve between the disordered phase and the ordered- 
uniform phase for ^ = 16 for 6i = —0.64. 
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Figure 6.6: Transition curve between the disordered phase and the ordered- 
uniform phase for ^ = 16 and the optimal exponent 6i = 0.55. 

We have evidence to believe the exponent 6i depends on R. But for the moment 
we content ourselves with choosing 6i as 0.64 and we consider the error as 0.30, 
this is 6i = 0.64 ± 0.3. 

Now that we have the collapse on N we can study the coefficient fi{R) in 
eq. K3). 



Step 2 
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Figure 6.7: Coefficients fi{R) for R = 2,A, 8, 16, 32 and the fit in eq. Q 
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6.1. Phase transition disordered to ordered-uniform 



In the previous figure 6.7 we estimate the errors on the fit of fi{R): 



f^(R) = (-0.31 ± 0.1)i?-(°-'^^^°-i). 



(6.7) 



Then the equation for the transition curve from the disordered phase to the 
ordered-uniform phase is: 



jY0.64±0.3 

fh^ = (-0.31 ± 0.1) A. 



(6., 



^0.64±0.1 

In terms of the constants A, B,C and D in eqs. (4.8)-(4.11) we re-writ^eq. (6.8) as 

C 



B ~ -2.01 



(AD) 



0.32 • 



(6.9) 



Checking the collapse on R for = 12 we obtain figure 6.8 
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Figure 6.8: Transition curve from the disordered to ordered-uniform phase 
for N = 12. 



We observed for R = 64 - where ~ = 0.1875 <C 1 ~ that the collapse given by 



eq. (6.10) is not good. 



Now we can check the collapse on R for N = 23, for ^ > 1, see figure 



6.9 
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Figure 6.9: Transition curve from the disordered to ordered- uniform phase 
for N = 23. 



Finally the collapse of data in eq. (6.10) considering all our data is shown in 
figure [67TO 
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Figure 6.10: CoUapse of data for the disorder to the order-uniform phase 
transition. 



We observe that for ~ < 0.375 and ^ > 8 the collapse is not goocj^ 

^The case N ~ 23, R — 4 could fall in this class, but is not possible to conclude due to 
the poor resolution for these points. 
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6.1. Phase transition disordered to ordered-uniform 



We conclude that the range where eq. (6.10) gives a good approximation is 

we have 

N 



f G (0.375 



If we want to include transition points for ^ ^ (0.375, 



to reconsider fi{N,R) defined in eq. (6.3) as a more complicated function on 
For the moment we exclude the data that are not in the interval (0.375, 8) and we 



present the figure 6.11 
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Figure 6.11: Transition curve from the disordered phase to the ordered- 
uniform phase. 



To conclude the present section, the expression of the collapse in figure 6.11 



IS 



ml = (-0.31 ±0.1) 



jY0.64±0.3 
^0.64±0. 



-A. 



(6.10) 
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6.2 Phase transition disordered to ordered non- 
uniform 

First we explore this transition curve for = 12, ^ = 8. 
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Figure 6.12: Transition curve from the disordered phase to the ordered non- 
uniform phase for = 12, ^ = 8. 



We observe that the transition curve shows curvature. The most natural fit we 
can propose is a polynomial where the coefficients are functions that depend on R 



and A^. In figure 6.12 we used a polynomial fit of degree 4: 



y{N = 12, R 



-0.09025 - 0.2389A + 0.0045A^ 



6.19 X lO'^A^ + 3.12 X lO^'A 



7\4 



(6.11) 



In the figure 6.12 we covered a large range of A. Nevertheless, to predict the 
triple point we can concentrate on "small" values of A but above Xtripie- For this 
range of values the transition curve can be approximated by a polynomial of a 



smaller degree. In figure 6.13 we present a linear fit for an interval of figure 
1632]: 



56 



6.2. Phase transition disordered to ordered non-uniform 
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i-m ■ ^ 
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Figure 6.13: An interval of the transition curve from the disordered phase to 
the ordered non-uniform phase for iV = 12, ^ = 8. 



First we explore a linear fit for the transition curves. In the subsequent section 6.4 
we will compare it with the results obtained for a fit using a polynomial of second 
degree. 



The fit of figure 6.13 



is: 



fh^ = -(0.11 ± 0.02) - (0.22 ± 0.01)A. 



(6.12) 



For general R we propose a linear fit for the transition curve from the 
disordered phase to the ordered non-uniform phase of the form 



fr? = ho{N, R) + hi{N, R)X 



(6.13) 



and we further make a factorisation ansatz and hi{N, R) = hi{N)hi{R), i = 0,1. 

First we concentrate on the scaling in N, keeping R = 8 and varying N = 
12,16,23,33. 
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NX 

Figure 6.14: Transition curves from the disordered phase to the ordered non- 
uniform phase for R = 8 



The fit for figure 6.14 is: 



-0.0178 ± 0.0022iVA - 0.125. 



(6.14) 



Here we conclude hi{N) cx N, ho{N) = const. 

Now we want the stabilise the triple point, at least the value of Xt a^we did 
in the previous section. Then if we re-scale the x-axis in figure 6.14| by a factor 
jYO.64 gg|. i-j^g figure 6.15 
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Figure 6.15: Approximate collapse of the slope for the transition curves from 
the disordered phase to the ordered non-uniform phase for R = 8. 



^To stabilise the value of m'^T it is not sufficient to re-scale the y-axis. 
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6.2. Phase transition disordered to ordered non-uniform 



We fit the transition curves for R = 4,8,32 for different by eq. (6.13) in 
order to find the coefficient hi{R). 
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Figure 6.16: Approximate collapse of the slope for the transition curves from 
the disordered phase to the ordered non- uniform phase for R = 2. 
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Figure 6.17: Approximate collapse of the slope for the transition curves from 
the disordered phase to the ordered non-uniform phase for R = 4. 
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Figure 6.18: Collapse of the slope for the transition curves from the disor- 
dered phase to the ordered non-uniform phase for R = 16. 
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Figure 6.19: Collapse of the slope for the transition curves from the disor- 
dered phase to the ordered non-uniform phase for ^ = 32. 

Now that we get the collapse on as hi{N) ~ we can determine the 
coefficients, we propose: 

hi{N,R) = -Nhi{R) (6.15) 



Figure 6.20 shows hi{R) and its fit. 
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Figure 6.20: Coefficients hi{R) for = 2, 4, 8, 16, 32 and the ffi OMAR-^-^'^. 
Then we have 

hi{N, R) = -(0.064 ± 0M7)NR-°-^^^^-^. (6.16) 



rr)2 _|_ 12.7 



Eq. (6.16) is the slope of the coexistence hne from the disordered phase to the 
ordered non-uniform phase. To solve the equation of this coexistence line we use 
the form: 

{m'^-m'^T)=hi{N,R){~X-~XT). (6.17) 

If we want to collapse the transition line completely, it is not enough to re-scale the 
y-axis. It is necessary to shift m^, for example, by substituting m-^ 
in figure 6.15 This leads to figure |6.21[ 
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Figure 6.21: Collapse of transition curve from the disordered phase to the 
ordered non-uniform phase for R = 8. 
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As another example for R 
-0.05 



4, form figure 6.17 we obtain figure 6.22 
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Figure 6.22: Collapse of transition curve from the disordered phase to the 
ordered non- uniform phase for R = 4. 



We conclude that for A around Xt the transition curve from the disordered 



phase to the ordered non- uniform phase obeys eq. (6.18) 

,2 _ ncA ^ n m^\i\Tf>~0.64±0.1 



-(0.064 ± 0.017)7Vi?^ 



A+ (2.697V°-^^ - 12.7) R 



1.92 



(6.18) 



6.3 Stabilising the triple point 

As we mentioned in section 16.11 we have the freedom to re-scale both axes of the 
phase diagram by a common factor. We want to use this factor to fix the triple 
point, i.e. we want to find a function of N and R, k{N, R), such that A;(A^, R)\t = 
const. 

To identify the triple point (At, {R, N) our strategy is the following: 

1. First we extract the dependence on R varying R and keeping fixed. What 
we would expect is of the form 

At = Z{N)R^^^^l (6.19) 
= M{N)R^^^^\ (6.20) 

For each N we estimate the exponents ei(A),e2(A) and the coefficients 
Z{N) , M{N) using the package "gnu-plot" . 
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2. In general the exponents can depend on N, but we wih show they do not 
and fix ei, 62 to a certain value. Since our final expression of the triple point 
may strongly depend on the choice of ei, 62, we will compare the expressions 
for two different sets of exponents]^ Both fits for ei in the case of - or 
62 in the case of m^j- - will be presented in the corresponding figures. We 
include the case when we vary ei and Z{N) - or 62 and M{N) in the case 



of m?" 



3. We extract the dependence on N proposing 

Z(N) = const. ■N'^\ 
M{N) = const. ■N'^\ 



(6.21) 
(6.22) 



Step 



First we explore the dependence on R. To do this we consider = 12 and 
several values for R and we estimate the intersection of both transition curves. 



Figure 6.23 is an example of this procedure. 
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Figure 6.23: Transition curves for iV = 16, i? = 8 



•X. 



The linear fit for figure 6.23 is given by the equation: 

ff? = -(0.090 ± 0.0003) - (0.0191 ± 0.0013)A. 



(6.23) 



^We will choose ei — —1.28, 62 = —1.92 vs. ei — —1.25, 62 = —1.89 and its correspond- 
ing coefficients Z{N),M{N). 
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The intersection of both curves is estimated in \t = 0.528ib0.02, 7n?T = — 0.2494ib 
0.003. 

Step 1 



Figure 6.24 shows the estimated values of for different R fixing = 12 



and a fit for these points: 
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Figure 6.24: Estimation of At for = 12 and three different fits. The fits 
were obtained via a function A^ = Z{N = 12)R^^ . In the first case we fit both 
parameters Z{N = 12), d and we got At = {Q.2Q±0M)R^-^-^^^°-°^^l In the 
second and third case we fixed the exponent Ci and vary Z{N = 12). For 



the second fit we chose ei = —1.28 and we fit Z(N = 12) as in eq. (6.24) and 
finally in the third case we fixed d = -1.25 to get Z{N = 12) = 8.34 ± 0.27. 

The fit for Figure [6?24l is: 

At = (8.20 ± 0.31)i?"^-28^°-\ (6.24) 

Then, if we multiply both axes in the phase diagram by a factor on R such that 
the X-axis is we stabilise the triple point to the value of At = 8.20(31). 



We can check the consistency of m^T with the eq. (6.10): 
^ We tried another exponent and we got an acceptable fit for At = 9.^1R^^^^^ . 
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Figure 6.25: Estimation of for = 12 and three different fits. The fits 
were obtained via a function Xt = M{N = 12)R'^'^. In the first case we fit 
both parameters M{N = 12), es and we got m| = (-10.24±0.66)i?-^-^^^°-°^. 



In the second case we fixed 62 = —1.92 and we fit Z{N = 12) as in eq. (6.25) 
and finally we fixed 63 = -1.89 to get Z{N = 12) = -13.06 ± 0.33. 

The equation of the cohapse is: 

= (-12.61 ± 0.20)i?(-^-92±o.i)_ (g_25) 

Then, if we re-scale the y-axis by a factor R^'^'^ we stabihze the value of R^'^^m'^T = 
-12.61. 

We conclude that for = 12 the equation of the triple point reads 

AT(iV = 12) = (8.20±0.31)i?"^-2^=^°-\ (6.26) 
mf.(7V = 12) = (-12.61 ±0.2)i?("^-92=^°-^). ^q^j^ 



Now, to estimate the dependence on we apply the same procedure as in the 
case N = 12, for = 8, 16 and A^ = 23. As in figures 6.24 6.25 we will compare 
three different fits. 
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Figure 6.26: Estimation of Ay for = 8 and three different fits. The first fit 



isAT(A^ = 8) = (9.20±0.9f)i?-^-^^^°-°'^^ The second fit is given by eq. (6.28). 
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Figure 6.27: Estimation of m^^ for N = 8 and three different fits. The first 
fit is m^T{N = 8) = -(f 0.9f ± 0.88)R"^-'^^^^-^^^. The second fit is given by 



eq. (6.29) 



We conclude that for N = 8 the equations of the triple point reads 
~Xt{N = 8) = (12.03 ±3)i?-i-28±o.i7^ 
fnUN = 8) = (-12.91 ±2)i?-i-92±0-22. 



(6.28) 
(6.29) 
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Figure 6.28: Estimation of At for = 16. The first fit is Xt{N = 16) 
(2.87 ± 0.60)i?-°-^^^°-^25_ rpj^g second fit is given by eq. (|6.30[). 
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Figure 6.29: Estimation of m^^ for N = 16. The first fit is rn^fiN = 16) 
-(4.9263 ± 0.91)i?-i-^^±°-^2^ The second fit is given by eq. (|6.3l[). 



At = 16 the equations of the triple point are 

At(A = 16) = (7.71±5)i?-i-28±o.43^ 
m^(A = 16) = -(9.76±5)i?-i-92±o.48_ 



(6.30) 
(6.31) 
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Figure 6.30: Estimation of for = 23. The first fit is Xt{N = 23) 
(10.77 ± 0.39)i?-i-5°^° °24. The second fit is given by eq. (|6.32[). 
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Figure 6.31: Estimation of m?T for = 23. The first fit is m'^T{N = 23) = 
-(18.73 ± 0.055)i?-2. The second fit is given by eq. (|6.33|). 
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At = 23 the equations of the triple point are 



Xt{N = 23) 
mj^{N = 23) 



(7.68 ± 3)i^"^•2'^=^°■2^ 
-(16.69 ±2)i?-i-92±o i2. 



(6.32) 
(6.33) 



We summarise the coefficients Z{N) obtained in eq. ( 6.26 )-( 6.33) in Table 6.1 



Table 6.1: Coefficients Z{N) and M{N) considering 



N Z{N) 



M{N) 



8 12.03 ±3 

12 8.20 ±2 

16 7.71 ±5 

23 7.68 ±3 



-12.91 ±2 
-12.61 ±2 
-12.42 ±5 
-12.92 ±4 



Figure 6.32 shows the values of Z{N) in Table 6.1 and three different fits: 
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Figure 6.32: Estimation of Z{N) for N = 8, 12, 16, 23. 



For m^T we present a constant fit in figure 6.33 
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Figure 6.33: Estimation of M{N) for N = 8, 12, 16, 23. 



Then we get as estimation for the triple point the eqs.: 

At = (45.45 ± 15)iV-0-64±o.2i;^-i.28±o.2^ (g_34) 
= -(12.7±l)i?-i-92±0-2. (6.35) 



In eq. (6.35) we observe that m^r seems not to depend on N. 

Let us consider the expression for the triple point fixing the exponents to 
ei = -1.25,62 = -1.89. The coefficients Z(iV)and M{N) are slightly different to 
those obtained with ei = —1.28,62 = —1.92, they are summarise in the Table 



Table 6.2: Coefficients Z(A^)and M{N) considering {\T{N),m'^T{N)) = 
(Z(iV)/?-i-25,M(Ar)i?-i-89). 



N Z{N) M{N) 

8 11.0 ±3 -12.91 ±2 

12 8.34 ±2 -13.06 ±2 

16 6.82 ±5 -12.42 ±3 

23 6.45 ±3 -12.92 ±4 



70 



6.4. Testing the fit of tlie transition curve I — II 



The values of Z{N) in Table 6.2 can be fitted via the function Z{N) 
(37.39 ± 30)iV"°-59=^°-2i, ^j^gj^ ^j^g expression for At reads 



At = (37.39 ± 3o)iV-0-59±0-2iii-i-25±o.2_ 



(6.36) 



For m?T we observe the values of M{N) in Table 6.2 are practically the same 
that in Table 



6.1 



We take as prediction of m?T eq. (6.37): 
m2T = -(12.7±l)ii-i-89±0-2. 



(6.37) 



For simplicity from all the fits presented in this section we choose those where 
the exponents can be written as integer multiples of 0.64. 



{Xt, ml) = ((41.91 ± i5)iv-0-64±o.20^-i.28±o.25^ _(i2.7 ± l)i?-i-92±0-20) 



(6.38) 



Note eq. (6.38) is in agreement with eq. (6.10) since if we substitute At in 
eq. (6.36) into eq. (6.10), the expected fh\ turns out to be: 



12.99i?^^-^^ 

2 



(6.39) 



eq. (6.37) coincides with the expression for in eq. (6.38) within the errors. 



6.4 Testing the fit of the transition curve / — 
// 

In this section we want to check the viability of a linear fit for the disordered to 
ordered non- uniform phase transition curve for values of A slightly above At. 
We study the significance of a quadratic term in A. We propose 

m\ = go{N, R) + gi{R, iV)A + g2{R, N)~X^ (6.40) 

From the previous section we adopt the assumption go{N,R) oc go{R) and 
gi{N,R) oc Ngi{R) and we propose g2{N,R) = g2{N)g2{R). 

We observe g2{N) = const. To find gi(R),i = 0,1,2, we keep = 12 fixed 

is an 



6.23 



and fit the curve (6.40) for the set of data with different R. Figure 
example for the procedure. 

Now we present the plots and fits for ho(R), hi{R) and gi(R),i = 0, 1, 2. 
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Figure 6.34: Coefficients /lo(-R) defined in eq. (6.13) and go{R) in eq. (6.40) 
for R = 4,8,32,64, 100. We can observe ho{R) ^ goiR) for R > A. The fit 



for go{R) is given in eq. (6.41). 



The fit for figure 6.34 is : 



go{R) = 18.33i? 



-2.63 
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Figure 6.35: Coefficients hi{R) defined in eq. (6.13) and gi{R) in eq. (6.40) 
for i? = 4, 8, 32, 64, 100. We can observe that for R > 4, hi{R) ^ gi{R). The 



fit for go{R) is given in eq. (6.42) 



The fit for figure 6.35 is : 



~gi{R) 



(6.42) 
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6.5. Collapse of observables 



As we can observe the coefficients obtained by a fit with a polynomial of degree 
1 and 2 are very similar. Thus the linear fit should be a good approximation for 
the region around the triple point 
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Figure 6.36: Coefficients g2{R) for = 4, 8, 32, 64, 100. 



(6.43) 



The fit for figure 6.36 is : 

g2{R) = 0.0014i?-2-29 + 1.97607 • IQ-^i?^-^^. 
The final expression for this transition curve is: 



m2 = (18.33i?"2-^^ - 10.2576^"^-^°^^^) - {0.195R~^-^^ + 0.004i?-°-°^^)iVA 



+144 • (0.0014i?-2.29 _^ 1.97607 x lO^'^ R^-'^'^)\^ + 0(A^ 



(6.44) 



6.5 Collapse of observables 



6.5.1 Collapse for A < Ay 

For this range of A we find that rescaling the x-axis by the factor A^^ the suscep- 



tibilities xo and Xi collapse. We show an example in figure 6.37 



^We chose the fit of the form ho{R) — aiR^^ + a2R^^ since we could not find a fit of 
the form ho{R) — aiR'^^ good enough to reproduce our set of data. 
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Figure 6.37: Collapse of xo and xi for R 
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Figure 6.38: Collapse of {(t>lu) for R = 8, N = 8. 
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Figure 6.39: Collapse of the Specific Heat in eq. (4.35) for R = 8, N 



We can see from figures 6.37 6.39 the equivallence between the Specific Heat 
criteria and the susceptibilities criteria at this regime. This subject is more widely 
disccussed in appendix [Cj 
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6.5. Collapse of observables 



6.5.2 Collapse for A 

If A is sufficiently large we expect that the relevant contributions to the action in 



eq. (4.12) are those of the potential terms. Depending on the value of R one or 



both of the kinetic terms (fuzzy and time derivatives) can be negligible. 

This leads us to reduced models, i.e. models that effectively depend on less 



parameters than those in our model in eq. (4.12) 



If the time derivative terms are negligible we would have a chain of A^-indepen- 
dent fuzzy spheres. If the fuzzy kinetic term is negligible the reduced model is a 
chain of matrix models, see Refs. |36]-|57]. If both kinetic terms are negligible two 
parameters are redundant, this lead us to the 1-matrix model (see Ref. |49] ) . 

We devote this section to the study of our model for A » Ay. We concentrated 
our analysis in the collapse of the phase transition and some observables. 

The transition curve for this region of the space of parameters can be well 



fitted by the expression (6.45): 



ml = wo{N, i?)A" 



(6.45) 



The exponent wi typically oscillates in the range [0.5, 1] as it can be appreciated 
in Table D.l in appendix [D| 

As an example we present the transition curve disordered — ordered non- 



uniform for N = 8, R = 16 and its proposed fit in Figure 6.40 
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Figure 6.40: Transition curve from the disordered phase to the ordered non- 
uniform phase for = 8, ^ = 16. The value of the exponent wi in eq. (6.45 ) 
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A characteristical behaviour of the specific heat (4.35) for the this region of 



large A is shown in figure 6.41 for A = 625, the last point in figure 6.40 (A^A 
5000): 
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Figure 6.41: Specific heat at A = 625, i? = 16, = 8. The critical point is 
estimated at = — 9 ± 0.5. 



The corresponding internal energy (4.34) and its partial contributions (4.40) 



(4.43) to figure 6.41 are shown in figure 6.42 
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Figure 6.42: Internal energy E in eq. (4.34) and its partial contributions 
eqs. (|I40|-(|I43| at A = 625, i? = 16, A^ = 8. 
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6.5. Collapse of observables 



Figure 6.42 shows that the leading contributions to the internal energy is 
from the potential. Although small, the contributions of the kinetic terms are not 
negligible. 

The next step would be to try to collapse in N and R the transition curves in 
this region of large A. As we have mentioned, the exponent wi typically oscillates 
in the range [0.5, 1]. The optimal value of wi strongly depends on the considered 
region of A -see Table 



D.l 



in appendix |D[ 

In section 6.2 we focused in the region around the triple point, there we con- 
sidered the transition curve a straight line, i.e. wi = 1. As we increase the range 
of A the exponent wi = 1 decreases. 



The iV-matrix model 



If the fuzzy kinetic term is negligible we arrive at a chain of matrix model. 

Ref. [47! studied the large N limit of this model, with a potential g(f)^. There 



a phase transition was predicted at the critical value (6.46) 



(-;.2)3/2 
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(6.46) 



where g is the critical value of the coupling and is the squared mass parameter. 



Under the appropiate translation to our parameters in eqs. (4.44)-(4.46), eq. (6.46) 
reads 



(6.47) 




In terms of the coefficient wo{N, R) and exponent wi in eq. (6.45 ) the prediction 



for the transition curve disordered ordered non-uniform at large A is: 



wo{N,R) 



/3iV2x2/3 2 



(6.48) 



The transition in figure 6.40 obeys eq. (6.47) (the fit in figure 
-0.12lp/3 and the prediction from eq. (|6.47|) is ml = -0.130A2/^). 



6.40 



IS m: 



Some other examples that have the transition in eq. (6.47 ) are shown in figures 
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Figure 6.43: Transition curve from the disordered phase to the ordered non- 
uniform phase for A^ = 16,^=16. 
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Figure 6.44: Transition curve from the disordered phase to the ordered non- 
uniform phase forA^ = 23,-R=16. 



In figure 6.44 we observe that the fit = — 0.45A^/'^ works well for A > 20 (the 



prediction for ti;o(23, 16) in eq. (6.48) is —0.53). Therefore if we only consider (or 



measure) A > 20 for this case the estimated value for wi would he wi > i. This 



can explain the different estimated values for wi in Table D.l in appendix [D| A 

where < | is that we 



D.l 



second possibility is that for those cases in Table 
have a transition of a different nature. This is analysed in the following section 
16X21 
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6.5. Collapse of observables 



The 1-matrix model 

In the 1-matrix model a transition is expected at 



-2VNX, 



(6.49) 



according to the notation in Ref. [30j where it was confirmed numerically, or 



-2VNc 



(6.50) 



according to the notation in Ref. [32] -see eq. (3.21 ). 



The question is if our model has such a transition. Theoretically for very large 
values of A and small value of R, the dominant term in the action is the potential. 
Then the model should effectively depend on less parameters. Following Ref. [50] 
we obtain the relevant parameters in our model form those in the 2-dimensional 
model: 



rrio 



'■2d - ^"T-Sd' ^2d = NXsd- 
As a next step we verify if our model has a transition at 



fhi (xV NX. 



(6.51) 



(6.52) 



With this intention we present figures 6.45 6.46 
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Figure 6.45: Phase transition to the disordered phase for several values of 
and R = 8. 
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Figure 6.46: Phase transition to the disordered phase for several values of 
and R = 16. 



For the points with larger A in figure 6.45 the = 16 data can be fitted by 

) = -0.52\/M + 5.8, (6.53) 



mi{R 



a line eq. (6.531: 



but a substantial difference is that it does not cross the origin as in eq. (6.52 1 
A similar situation occurs for ^ = 16 in figure 



6.46 



for A^ = 16, 23 where the 



fit is given by eq. (6.54): 



fhl{R = 16) = -0.26 VA^A + 2.8. 



(6.54) 



Note that the coefficients in eq. (6.53) for R = 8 are approximately doubled com- 



pared to eq. (6.54) for R = 16. 



We conclude that we cannot confirm the phase transition in eq. (6.51) for the 



2-dimensional model. A possible explanation is that in the data obtained we have 
not reached a sufficiently large A. 

Now we focus on the collapse of other observables. 



Collapse of 0^^^ 

In this section we investigate the collapse of the norm of the field in eq. (4.27) for 
large A. This quantity is of interest since for the ordered non-uniform phase for 



80 



6.5. Collapse of observables 



large A the main contribution to {4>iu) comes form higher modes, and the other 



quantities we are measuring are related to the lowest modes. In addition it is 
known from the 1-matrix model studied in Ref. [30j that for a region where the 
kinetic term is negligible, (0^//) has the simple form: 



27rr 



(6.55) 



where r is the squared mass parameter. It would be interesting to check if our 
model, under the appropriate translation of parameters, can reproduce the 1- 
matrix result as a limiting case. 



For m? < m?c we found that the norm of the field, {(fP' 



all! 



eq. (4.27) 



be well fitted by a line. A important difference from the 1-matrix model studied 



in Ref. [30] is that the line does not crosses the origin. We propose eq. (6.56) 



t>lu) 



vo{N, R,X) + vi{N,R,X)m^. 



(6.56) 



A concrete example is presented in figure 6.47 for the same parameters as in 
figures [eTIT] - [61421 

0.014 



0.012 - 

0.01 - 

0.008 - 

0.006 - 

0.004 - 

0.002 - 
0^ 



-10 



{"Pall, 

-O.OOSlm^ - 0.0185 



-\ — I 
- X- - -; ' 



+ + 



¥ )K ^ )K • 

X X X X )( : 



-9.5 



8.5 

"2 



-7.5 



m 



Figure 6.47: (flu) and a estimation vs. at A = 625, -R = 16, = 8. 
We observe that the contributions from the zero and the first mode, {(Pq) 
and {(fil) respectively to (fan) are small as we expected. The critical value 



IS 



-9.0 ±0.5. 
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Note: There are some technical difficulties in the measurement of vi and vq: 
The coefficients vi and vq must be measured on an appropriate range of , where 
vi and Vq stabilise. This typically happens for < m'^c, i-e. in the ordered 
non-uniform phase where we have thermalisation problems. For the parameters in 



figure 6.47 they appear for m? < —15. Then in figure 6.47 the value of vi and 
Vq stabilise for —15 < < —8.3. 

In addition to the 1-matrix model, our model has another limiting case: the 
chain of matrix models. In the chain of matrix models the fuzzy kinetic term 
should be negligible while in the 1-matrix model both kinetic terms are negligible. 

We conjecture that the non-vanishing coefficient vo{N, R,X) is related to the 
chain of matrix models and therefore should reduce to zero for the limiting case of 
the 1-matrix model. Therefore the coefficients vq(N, R,X) and vi{N, R,X) should 
depend differently on N and R. 



Finally we present the attempts to collapse {(p'^.ji) for A > At in figures 
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Figure 6.48: Collapse of {(f)lii) ioi R = 8, N = 8. We re-scale the x-axis by 
the factor A"°-^. For values of A slightly above At the collapse works but not 
for A > 1500. 



In figure 



6.48 



the X-axis is re-scaled by the factor A''. The optimal value of tj 
depends on the range of A. rj = —0.8 is optimal for 40 > A > At, it gives and 
acceptable collapse for A < 100 in the range r] £ [-0.75,-0.85]. For 2500 > A > 
1500 we choose 77 = —1 as in figure 6.49 
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Figure 6.49: Collapse of {4>lii) for R = 8, N = 8. We re-scale the x-axis by 
the factor A^^. For values of A > 1500 the collapse is valid. 



We conclude that the collapse of (0^^^) depends on the range of A. For the 



Specific Heat we have a similar situation a for 
For A < 100 we have in figure 



6.50 



the collapse of the Specific Heat re-scaling 



the X axis by the factor A with rj = —0.8 as in figure 6.48 
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Figure 6.50: Collapse of the Specific Heat in eq. (4.35) at i? = 8, 
100 > A^A > 16. 



8 for 



In figure 
in figure 



6.51 



6.50 



we present the collapse of the Specific Heat for the same data that 
but re-scaling the x axis by the factor iV-2/3^4/3;^-2/3_ 
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Figure 6.51: Collapse of the Specific Heat in eq. (4.35) at i? = 8, = 8 for 
NX > 100. 
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Figure 6.52: Collapse of the Specific Heat in eq. (4.35) at i? = 16, iV = 8 for 
NX > 1000. 



From figures 6.51 6.52 we observe that re-scaling the x-axis by the factor 
^-2/3^^4/3^-2/3 £^ ^j^g critical value of to iV-2/3i?4/3^-2/3^2 p. -0.031. 

We have the same case in figure 6.53 for N = 12: 
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Figure 6.53: Collapse of the Specific Heat in eq. (4.35) at i? = 16, 
for A^A > 1000. 



12 



We conclude the collapse of observables and the transition curve for A > At 
does not lead to the 1-matrix model. Furthermore, the collapse of observables 
indicates that the exponent of collapse — rj for {if'^n) as in figures 6.48- 6.49 — 
depends on the range of A. 

For the largest values of A considered our results agree to those predicted for 
the chain of matrix models in Ref. [57] at large A^, as it is shown in figures 
[63T1I6331 



Chapter 7 

Discussion of the results 



As it was mentioned in section 4.3 we are interested in taking the limit — > oo. 
If we scale R in terms of N we can access different limiting models. The key point 
in our analysis is the behaviour of the triple point and the phase coexistence curves 
under those limits; it will decide which phases survive at the end. 

If the triple point remains finite, it indicates that the three phases exist. If the 
triple point goes to infinity it would indicate the existence of the Ising type phases 
while the ordered non-uniform phases would disappear. 

Now we proceed to analyse the model around the critical curves. Coming back 



to the discretised model given by eq. (4.12), it can be re- written as 



^* r 1 
5 [$] = ^ Tr I A$ (t) (t) + D + At) - $(t))2 + B^^{t) + C$^(t) (7.1 



t=i 



where we repeat for conveniece the definitions (4.8)-(4.11| 

27rAt 

N ' 

D 
B 
C 



NAt' 

N ' 
vri^^AAt 



N 

We have the freedom to re-scale the field $ in the following way: 



1 



(7.2) 
(7.3) 
(7.4) 
(7.5) 

(7.6) 
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and the constants (7.2) - (7.5) change to 

27rzAt 



A = (7.7) 

° = -NAT- <^-'*> 
B = (7.9) 

C = (7.10) 

As we mentioned in section 14.4.31 we can choose z such that we fix one of the 



constants (7.7)-(7.10). We fix yl = 27r, i.e 



N 
At' 



Under the re-scaling (7.11) the constants (7.7)-(7.10) change to 



A = 27r, 
D 
B 
C 



27ri?2, 



(Ai)2 
R 

■kR^XN 



At 



ttNR^X. 



(7.11) 

(7.12) 
(7.13) 
(7.14) 
(7.15) 



We re- write the obtained expression of the triple point of eq. (6.38) making 



the substitutions of eqs. (4.44)-(4.46) to get 



{AtXrdAtfm-) = (^41.91 {^)\ -^^-^ {^f^ 

with 7 = 0.64 ± 0.2. 

We scale i? = iV^ for /? > 0, then R — > cxD as — > oo. 



(7.16) 



Eq. (7.16) changes to eq. (7.17) 



(AtAr, (At)2m^) = (41M{At)'^^N-^'^^'^+^\-12.7{Atf^N- 



■37/3 



Choice of At. We choose 



At 



1 



(7.17) 



(7.18) 



with K > as an ansatz. 



7 Discussion of the results 



87 



Then the time extension T is 



NtAt = N 



l-K 



(7.19) 



for Nt = N. 
If K < 1 then T 



cxD as 



CXD. 



Implementing eq. (7.18) in eq. (7.17) we get 



(Ar,m|) = (AtAT,At2m2) = (41.91iV-^(2/3+2«+i)^ _^2.7iV-37(/3+«)^ . (7.20) 



We note that the dimensionless quantities in eq. (7.20) goes to zero as N — > 00 



Eq. (7.20) can be written as 



{XT,ml) = (41.917V'=(i-2^)-^(i+2/^),-12.7iV2K-37(/3+.) 



(7.21) 



To analyse the limit — > 00 we focus on the exponents of N in eq. (7.20). For 



Xt the exponent is k(1 — 27) — 7(1 + 2/3) and this is clearly negative if 1 — 27 < 0. 
Although the error in 7 is quite large even taking the lower bouncQthe exponent 
k(1 — 27) — 7(1 + 2/3) is negative, which implies 



0, N 



00. 



For the exponent of iV is 2k; — 3^{(3 + k), it is positive or negative depend- 
ing on the values of /? and k. 

In any case what we would have is the disappearance of the ordered uniform 
phase and therefore the phase diagram of the commutative theory cannot be re- 



covered from the one of our studied model (7.1 ) 
The tricritical action read^ 



St['^,N,R] 



N 2 

Y: Tr [^'^ it) (t) + ^ f 1,1 mt + At) - mf 

t=l 



.N 
25.47r 



N \AtJ 



N 



$'(t) + 



41.97r 



N \NAt 



(7.22) 



We consider the particular case At 



and i? = N^, 

/N ' 



1 The lower bound would be 0.64 - 0.2 = 0.44. 
^We chose 7 « |. 
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7.1. Comparison with other numerical studies 



N 



N 



(7.23) 



The exponent of for the $ -term in eq. (7.23), 7(/3 



1, is negative for 



the values of /? considered in the limits in section 4.3 -see eqs. (3.13)-(3.15 ). We 



observe that for N large the leading contribution to eq. (7.23) comes from the 
temporal kinetic termj^ 

If we compare the powers on N multiplying the kinetic terms in eq. (7.23), we 
see that the temporal kinetic term is A^^^^-*^ times larger than the fuzzy kinetic 
term. Under the limit N — > cxd the fuzzy kinetic term is negligible and then, 
because in some sense the geometry of the sphere is screened by the temporal 
kinetic term, the model behaves more like a "matrix chain" system interacting to 
first neighbours with the potential Xcj)^ (see Ref. |46j). 

To maintain the uniform ordered phase in the limit — > oo it is necessary 
to reinforce the fuzzy kinetic term as in Refs. l22]-[29]. We will come back to this 
point in the conclusions. 



7.1 Comparison with other numerical studies 



7.1.1 on the fuzzy sphere 



We found the 3 phases present in the model in eq. (3.20) -see Refs. |30j-|32j. 
A substantial difference from the 2-dimensional case studied in Refs. I30I-I32I 



is that the whole phase diagram in figure 3.1 collapses using the same scaling 
function of A^. We do not have this situation since the scaling on A^ for the 



the transition curves (6.10) (6.18) are different. The collapse in the 2 dimensional 



studies do not depend on the radius since it can be absorbed in the couplings. This 
is in contrast to the 3 dimensional case, where the radius plays an independent 
role. 



In section 6.5.2 we studied the collapse of observables and of the transition 
curve disordered to ordered non-uniform. We conclude that we do not observe 

3We will identify J2t Tr (t) (*)) / Et Tr [^(t + At) - <^{t)f. 
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the same scaling behavioui]^ that as in the 2-dimensional case. Furthermore, the 



collapse of the transition curve / — /// for large A - see eq. (6.47)- indicates that 



in this regime our model behaves as a chain of interacting matrices, where the 
fuzzy kinetic term is neglected Ref. |47) . 



7.1.2 Non- commutative lattice studies 

In this section we want to compare our results with those obtained in Refs. |35].[36l 
and j37| where the following model was studied: 



T 

E 

t=i 



E 



m 



D,^{t)D^ 
A 



t 



Ht)f + \m+i)-my 



2 ^\t) + -^\t) 



is an hermitian matrix. Eq. (7.24) describes a scalar field 
commutative torus and interacting under the potential. 



(7.24) 
living on a non- 



Comparing eq. (7.24) to eq. (7.1 ) we observe that, besides the different interpre- 



tation in each discretisation scheme, the substantial difference between the models 



is in the spatial kinetic term. The first term on the right-hand-side in eq. (7.24) 



corresponds to the energy due to spatial translations on a squared N x N lattice 



with lattice spacing a, while the first term on the right-hand-side in eq. (7.1) cor- 



responds to the energy due to translations (or rotations) over a discrete version of 
the sphere. It is to be expected that in the regime where the spatial kinetic term 
is negligible, both models describe the same physics. Now we want to compare 
the phase diagram for each model, but first we denote the constants of the model 



given by (7.24) as follows: 



A 
D 
B 
C 



N 
2"' 

2"' 

2 ' 
NX 



(7.25) 
(7.26) 
(7.27) 
(7.28) 



*This is under the appropriate translation of parameters from the 2-dimensional model 



to the 3-dimensional model, see eq. (6.51 1 
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7.1. Comparison with other numerical studies 



In Refs. [35] -[36] the phase diagram shows the existence of three phases: 

• Disordered phase 

• Uniform phase 

• Striped phase 

This phase diagram stabilises taking the axes as A^^A vs. N'^m?, this is, under 
the re-scaling of the axes A — > N'^X and — > N'^m? the transition lines collapse 
as follows: 

• Disordered phase — Uniform phase coexistence line: 

N'^rr? ^ -O.SSiV^A. (7.29) 

• Disordered phase — Striped phase coexistence line: 

iV^m^ ^ -0.52iV2A - 64, (7.30) 

and the triple point is given by 

(iV^A, N'^m^) ^ (220, -150). (7.31) 



In the model (7.1 ) we found the existence of three phases: 



I: Disordered phase 

II: Ordered uniform phase 

III: Ordered non- uniform phase 



We conjecture that the striped phase of the model in eq. ( 7.24 ) corresponds to the 



ordered non- uniform phase in the model given by eq. (7.1). 



Coming back to the fuzzy model in eq. (7.1), their coexistence curves I — II 
and / — /// stabilise under different scaling: 

• Disordered phase — Ordered uniform phase: 

i?37^2 = _o.31Af^i?27X. (7.32) 
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Disordered phase — Ordered non-uniform phase: 
12.7 



0M4{NR^yX + 2mN^-°-^^R^'^ (7.33) 



The triple point is given by the equations: 

{XT,m^T) = (41.91iV-0-6"±°-20i?-i-28±0-25,_(i2.7±l)i?-i-92±0-20) 



(7.34) 



In the region around the transition curve I — II (to the ordered uniform phase) 
the kinetic term is relevant. Then, because of the difference in the nature of the 



kinetic term in the models in eq. (7.24) and eq. (7.1), we cannot expect that there 



exists a re-scaling such that both collapses are compatible. In the region around 
the transition curve / — III (to the ordered non-uniform phase) the leading term 
is the one due to the potential, then we could expect that there exists a re-scaling 
such that both collapses are compatible for the transition curve / — III. Now we 
consider an special case of R and At. 



Case R = N, At = l. 



We consider R = N and At = 1. Implementing the substitution in eq. (7.32) we 
get 

-0.31iV^^A. (7.35) 



Comparing eq. (7.29) to eq. (7.35) and considering 37 « 2, we conclude that the 



transition line I — II for both models collapse with the same dependence in A^. 
We have the same situation when we compare the triple points. Implementing the 



substitutions R = N and At = 1 in eq. (|7.16|) we get 

1 



(AT,m^) = 41.91 



-12.7 



1 



Ar37 



41.91 12.7 



Ar2 ' 



(7.36) 



The dependence of the triple point on is the same as in Refs. [35] -[36]. The 
critical action is: 

N 



t=i 



-27r 
_iV 

25.47r 



[• 
(t) (t) + 27tN [$(t + At) - <^{t)Y 



N 



,9 , , 41.97r ,4, ■ 



(7.37) 



From (7.37) we observe that the dominance of the temporal kinetic term is even 



stronger than in eq. ( 7.23 ) 
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7.1. Comparison with other numerical studies 



• Observation: 



Choosing (3 = ^ in eq. (7.23) we get the same critical action as in eq. (7.37) 



but the dependence on in the triple point is not the same. 

1 is: 



The triple point ( |7.20[ ) for /? - 2 

(AT,m^) = f41.917V5-37^_i2.7Ari-37^ . 



(7.38) 



For the transition curve / — III the predictions from models ( 7.1 )-( 7.24) are 



different. The main difference is that the transition curve in the model given by 
(7.1) is a curve as it is shown in figures 6.43||6.44 whereas in the model given 



by eq. (7.24) this transition curve is a straight line. Nevertheless, we found that 



in a range of parameters around the triple point the transition curve / — /// in 



the model (7.1 ) could be approximated by a line as in (7.24). We conjecture that 



this is the range of parameters studied in Refs. [35]-[36]. 



Taking R = N and At = 1 in (|7.33|) we get 

12.7 



0.064iV3TA + 2.77V^^-^-^\ 



If we consider just the terms in and A in (7.39), we obtain N"^^ 



m vs. 



(7.39) 



We observed 37 ~ 2 while 4j — 1 = 1.67 ±0.8. We conclude that the dependence 



on in the eq. (7.39) could be the same as in eq. (7.30). 



Chapter 8 

Conclusions from part |I| 



We presented a numerical study of the Xcp'^ model on the 3-dimensional 
Euclidean space which was regularised by means of: 

— the fuzzy sphere Sp for the spatial coordinates 

— a conventional lattice with periodic boundary conditions for the time 
direction. 

The obtained model was 

t=i ^ ' 



where ^{t) G Matj-^, for t = 1, . . . , iVj. 



We found the phase diagram of the model according to the specific heat. 
Following this criterion we determined the critical values of A and for 
fixed N and R denoted by an index "c" , Ac and m^. In addition we found the 
tricritical point (tti^jAt)- At divides the phase diagrams into two regions 
according to the behaviour of the observables. 

— At > A > 0. In this domain we observe Ising type orderings. For 
vr? > ?Ti^ we have a disordered phase. For w? < we found a 
uniform ordering. 
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— A > At. For > we have a disordered phase. As in the two di- 
mensional model discussed in chapter [s] we found for m? < there is 
a non-uniform ordering. For <C we encountered thermalisation 
problems, therefore we cannot conclude if there exists a boundary for 
the non-uniform ordered phase as in the 2-dimensional case |32j . 

We compared our results with those obtained by other criteria (two point 
functions of the different modes) and we conclude both criteria are quali- 
tatively equivalent. We followed the specific heat criterion which is a more 
universal quantity: it captures the phase transition without taking into ac- 
count which is the dominant mode. 

We found the existence of three phases: 

— /: Disordered phase. 

— II: Ordered Uniform phase. 

— Ill: Ordered Non-Uniform phase. 

These three phases were also found in the 2-dimensional A(/)^ model on a 
fuzzy sphere studied in Refs. [30], [32] . 

The phase of non-uniform ordering is characterised by the dominance of 
several angular momenta for / > 0. In this phase the rotational invariance 
is broken. 

We get the transition curves of the model: 

— I — II: It turns out to be a line given by the equation: 

jY0.64±0.3 

(At) V' = (-0.31 ± 0.1)— ogj^AtA. (8.2) 

— I— III: We can observe that the transition curve shows curvature. The 
most natural fit we can propose is a polynomial where the coefficients 
are functions that depend on R, N and At. Nevertheless, to predict 
the triple point we can concentrate on values of A slightly above Ay. 
For this range of values the transition curve can be approximated by a 
line. We conclude that for A around At the transition curve from the 
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disordered phase to the ordered non-uniform phase obeys the following 
equation: 

-0.64±0.1 

\At J 



{Atfml = -(0.064 ± 0.017) (^) AtX + 



(2.69iV°-^i-12.7) f^j . (8.3) 

— We conjecture the existence of a transition curve II— III as in Ref. [32] , 
but due to thermalisation problems it was not possible to measure it. 

— The effective action has several minima, and the thermalisation prob- 
lems appear when it is not possible for the algorithm to tunnel between 
those minima. We sketched the main technical features of these ther- 
malisation problems. 

We obtained the equation for the triple point: 

'1.92±0.20\ 



(AtAT,(Ai)2m|) = Ul.917V-°-«4±°-20i?-i-2™5,_(i2.7±l) (^^^ 



(8.4) 



• Different limits of the fuzzy sphere can be obtained scaling i? as a function 
in and taking the limit N — > oo, 

— Commutative sphere: = const., N — > oo 

— Quantum Plane: R'^ cc N 

— Continuum flat limit. It requires: 

^ 1 > e > 0. (8.5) 



We analysed the behaviour of the scalar model in (8.1) under the different 
limiting spaces above. 

Our numerical results reveal that the triple point goes to zero under the 
limit N — > oo, this is valid for all cases considered in the previous point. 

In other words, in the limit N — > oo the non- uniform ordered phase domi- 
nates the phase diagram. This result is as consequence of the UV-IR mixing: 
integrating out high energy in the loop produces non-trivial effects at low 
external momenta I20I-I21I. 
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A perturbative analysis of the action (our eq. 4.4 of chapter |4] ) : 



+ '^<^Ht) + ^<^\t) , (8.6) 



is presented in Ref. |29]. It shows that the action given by eq. (8.6) does 
not reproduce the commutative continuum hmit. This result was obtained 
from a expansion to two loops. Our results show that the UV-IR mixing is 
presented for all values of A. 



Following Refs. [22] and |29|j the action (8.6) should be corrected by modi- 
fying the action. 

This modification in the action is equivalent to reinforcing the fuzzy kinetic 
term and it can be achieved adding a term <I>(£^)^<I>/(A^i?^) where A is a 
momentum cutoff. 

We compare our results to those obtained in a numerical study of the A</)^ 
potential on a non-commutative torus (see Ref. |35j-|36j). Our results are in 
agreement if the parameters R and At are scaled appropriately. 

The axes cannot be chosen consistently for all regimes to lead to a sta- 
ble phase diagram for large A^. This is a significant difference to the 2- 
dimensional formulation which has the "privileged" property of leading to a 
phase diagram. 

We compare our results at large coupling A to those of the chain of matrices 
at the large A- limit in Refs. |47j-|64j (where this model is known as the c = 1 
model in string theory). Our simulation results fully agree. The disordered 
to ordered non-uniform phase transition in this regime obeys the predicted 
transition in the model c = 1 in the large A-limit. 



^To be precise, in ref. ^29j discusses S"!. x M. 



Part II 
Scalar Field Theory on S 
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Chapter 9 

Prom fuzzy CP^ to a fuzzy S 



5"^ is an especially important example since it is the most natural replacement of 
in studies of Euclidean quantum field theory. Therefore our motivation to propose 
a fuzzy approximation to S"^, namely Sp. Nevertheless cannot be quantised in 
the strict sense since it is not a phase space. It is important to clarify how we can 
we obtain a matrix approximation to it. 



9.1 Sj^ in analogy to S'^ 



Having in mind the seminal example of fuzzy space, we search for five matrix 
coordinates, Xa, a = 1 ... 5, which fulfil a matrix equation of a 4-sphere in M^: 

5 

Y,XaXa = R^l. (9.1) 

a=l 

We propose 

Xa = -^Ta, (9.2) 

with Ta G Mati the Dirac matrices including 75. They obey the algebra: 

{r^.r^} = 25,a. (9.3) 



Eq. (9.1) follows directly from eqs. (9.2) and (9.3). 

The next step is to propose a sequence of matrices approximating S^. To 
achieve this we note that Ta (which are 4-dimensional matrices) give the represen- 
tation (2, 2) of Spin{5). We can therefore consider the irreducible representation 
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9.1. Sp in analogy to Sp 



obtained from the L fold symmetric tensor product of this representation i.e. the 
Spin{5) representation j). It will contain a set of five matrices: Ja, a = 1 ... 5 
which can be realised as the symmetrisation of L copies of the T matrices in the 
Spin{5) fundamental representation: 



Ja = Ta®!®-- -(8)1 +l(8)ra® •••01 + ••• + 1(8)1(8) •••®ra . (9.4) 




L— terms 

The subscript sym indicates that we are projecting onto the irreducible totally 



symmetrised representation. The dimension of the matrices defined in eq. (9.4) is 



6 

and Ja G Mat^i^. They satisfy the relation: 

JaJa = L{L + A)l. (9.6) 



Now we can generalise the definition of the matrix coordinates (9.2) to eq. (9.7): 



VHL + 4) 



The definition of the matrix coordinates (9.7) guarantees that the matrix equation 



of a 4-sphere (9.5) is fulfilled for matrices of dimension d^. In the limit L — > oo 



the matrix coordinates in eq. (9.7) commute. In this limit we recover the algebra 
of function of continuous S^, C°°(5'^) . 

However, a substantial difference to the fuzzy 2-sphere case is that for finite 
L the matrix coordinates Xa do not provide a complete basis for the algebra of 
functions, i.e. Xa cannot provide a basis for Matdj^- 

To clarify this point we analyse the lowest approximation, L = 1. We define 
dab proportional to the commutators of Ta- 

aab = Yi[Ta,Tb]. (9.8) 

If is a matrix representing a function on 5^, it will be of the form 



F = Fol + FaTa. 



(9.9) 
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However, a matrix product of two functions of the type in eq. (9.9) will involve 



non-zero coefficients of a^h, the matrices in (9.8) 



F' = F^1 + F^r, + F>,6. (9.10) 



The 10 coefficients F^^^ in eq. (9.10) have no corresponding counterparts in the 
expansion of functions on commutative S^. In the language of Statistical Physics 
these coefficients constitute a set of extra degrees of freedom. On one hand, if we 
exclude them the involved algebra is not associative. On the other hand, if we 
included them, the approximated space is not exactly 5^. 

The first option, advocated by Ramgoolam [13], is to project out such terms, 
in which case one is left with a non-associative algebra. This involves additional 
complications and does not seem particularly suited to numerical work. In addition 
the necessary projector must be constructed. We will return to this point in chapter 



11 where we will, in fact, give the projector. 

An alternative is to include arbitrary coefficients of aab (demanding an asso- 
ciative algebra) and attempt to suppress such coefficients of unwanted terms, by 
making their excitation improbable in the dynamics. In this approach our algebra 
will be a full matrix algebra and obviously associative. 

Including the extra degrees of freedom will lead us to work in a bigger space. 



This bigger space is CP , in section 9.2 will review its construction. 



9.2 Review of the construction of CP^ 

The fuzzy version of CP^-i denoted by the subindex "F", CP^'S is a matrix 
approximation to the continuous CP^^"*^. In this section we review the construction 
of CP| following Ref. pT^. 

A standard definition on CP^~^ is the space of all norm-1 vector in modulo 
the phase. For any unit vector \ip) we can define a rank-one projector, 

Viij) := (9.11) 

Then CP^~^ can be defined as the space of all rank-one projectors 



CpN-i ._ |p ^ MatN] = V,V^ = V, TrV = 1}. (9.12) 
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9.2. Review of the construction of CP^ 



To construct the set of global coordinates for CP^~ we need a set of N"^ 
hermitian matrices {l,t^}, /i = 1 . . . A^^ — 1. The set {t^} is a basis for the Lie 
algebra of SU{N), normalised as 

Tr (t^U) = 6f,^. (9.13) 

For our case of CP^ we start with Spin{6) = SU{4). Let {Jab}, A,B = 1 . . .Qhe 
the Spin{6) generators. They are equivalent to the generators of SU{A). As a 
basis for our algebra we take the set {A/is}, A,B = 1 ... 6: 

Aab := -^^AB, (9.14) 

their algebra is: 

^AB^CD = A.AB;CD^ + ^^e^ABCDEF^EF + (9.15) 

{6aC^BD + ^BD^AC - SbC^AD - ^AD^Bc) ■ (9.16) 



2^2 

Aab-cd is defined as a two indexes Kronecker symbol 



A-AB;CD = - {Sac^bd - Sad^bc) ■ (9-17) 



The algebra given by eq. (9.16) admits representations of dimension 
di = for L integer. For the lowest non-trivial level L = 1 we have 

the explicit form of the generators in appendix |Ej 

A projector V can be expanded in terms of the basis given by eq. ( 9.16| ): 

V = ^1 + UbAab, (9.18) 

(,AB are fifteen real coordinates. 

= V implies that £^ab obeys the restrictions: 

UbUb = \, (9.19) 

^ABCDEPiABicD = S,EF- (9.20) 



Taking contractions of eqs. (9.19) and (9.20) we get the three identities: 
1 3 

-eABCDEF£,AB^CD£,EF = 7, (9.21) 



2V2 ' ' ' 4' 

^^^ABCDEF^AB = '^CabCcD - "^^AC^BD + "^^AdCbC, (9.22) 

UcCcB = -l^AB. (9.23) 
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Eqs. (9.19)-(9.20) describe how CF% is embedded on M}^. The global coordinates 
{Cab} allow us to describe the geometry of CP p. Following Ref. [18j we have the 
geometrical structures 



^AB;CD 
PaB;CD 

^AB,CD 
JaB;CD 



[PaB;CD ± iJAB;CD) , 



■^Aab;CD + v^c^fecD'^SF - '^CabCcd, 
-Aab;CD - V2d^^(jjj^EF + '^S.abCcd, 

^fABCD^EF- 



(9.24) 

(9.25) 

(9.26) 
(9.27) 



Kab;CD is the Kdhler structure, Jab;CD is the complex structure and Pab;C'D is 
the metric. 



From eq. (9.16) we get the explicitly form of the normalisation constants, 



eq. (9.28), and the structure constants, eq. (9.29): 



dABCDEF 
fABCDEF 



^ABCDEF, 



{SacAbD;EF — SadAbC;EF 

+SbdAac:EF — SbcAad:EF 



(9.28) 



(9.29) 



From here to the end of the thesis, we used the Kronecker's delta to arise and low 



indexes unless the opposite is indicated. A simplification for eqs. ( 9.25 )-( 9.27) is 
the following: 



Pab;CD 

^AB;CD 
JaB:CD 



^Aab;CD - 2 {iAcisD - CadCbc) , 

-AaB;CD + {£,Ac£,BD - iAoiBc) , 



V2 



{^acCbd - SadCbc + SbdCac - ^bcCad) 



(9.30) 
(9.31) 
(9.32) 
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9.2. Review of the construction of CP^ 



Some properties of the metric Pab;CD and the complex structure 

JaB:CD- 



PaB;CD 


= PcD;AB = —PbA;CD, 


(9.33) 




— ^CD;AB — ~^BA]CD^ 


(9.34) 


JaB;CD 


= —JCD;AB = —JbA;CD, 


(9.35) 


^AB;CD 


'■= PaB;EfPeF;CD = PaB;CD, 


(9.36) 




'■= PaB;EfPeF;CD = PaB;CD, 


(9.37) 


'^AB;CD 


'■= JaB;EfJeF;CD = —PaB;CD, 


(9.38) 


JaB;EfPeF;CD 


= PaB;EfJeF;CD = JaB;CD- 


(9.39) 



Pab;CD, Pab cd ^^'^ Kab\CD are projectors, their ranks are given in eqs. (9.40 ) 



(9.421: 



PaB;AB 


= 6, 


(9.40) 


P-L 

^AB]AB 


= 9, 


(9.41) 


^AB;AB 


= 3. 


(9.42) 



Pab;CD projects on to the tangent space of CP^ and P^b ab onto the orthogonal 
compliment in M}^. 

9.2.1 CP^ as orbit under Spin{6). 

Now we want to perform an explicit construction for CP'p, we will analyse the 
induced line element. CP'p can be obtained taking one fiducial projector and 
rotating it with and element of Spin(6): 

P{il;) = U{ip)P^U-^{ip), U{^) G Spin{6). (9.43) 

We choose P^ as 

P' = \l + eAB^AB 

= |l + -^(Ai2 + A34+A56). (9.44) 
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Now we are interested on calculating the line element, ds'^, at the north pole, 



defined by eq. (9.44). 



(9.45) 



AB 



Spin{6) rotates Aab as a tensor, 



£,AB = RacRbd£,cd- 



(9.46) 



The sum d^^^ in eq. (9.45) can be written in terms of an infinitesimal rotation 
R^^dR: 



Y,deAB = -Tr [R-'dR,e 

A,B 

where ^ is the matrix of coefficients and 



(9.47) 



R ^dR := -cabLab- 



(9.48) 



Eq. (9.48) is known as the Maurer-Cartan forms for the group of rotations on 



6 dimensions. Lab are the generators of this representation. 
The line element turns out to be: 

ds'^ = ^ [(ei3 - 624)^ + (ei4 + 623)^ + (ei5 - 625)^ + (eie + £25)^ 

+ (e35 - 646)^ + (636 + 645)^] • 



(9.49) 



From eq. (9.49) we corroborate that CP^ is a 6-dimensional space 



Tangent forms 

We define the tangent forms as 



Sab = PAB;CDeCD- 



(9.50) 



At the north pole in eq. (9.44) we have the non- vanishing coordinates are ^ 



12 



^;^,^34 = = and permutation of them. The tangent forms at the 

north pole turn out to be: 
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9.2. Review of the construction of CP^ 



-12 



'13 



'14 



'15 



'16 



-34 



-35 



-36 



-56 



(ei3 - 624) 

(ei4 + 623) = e 

(ei5 - 626) 

(ei6 + 625) = e 



-24' 



23' 



-26' 



25' 



(635 - 646) = -e. 



46' 



(636 + 645) = e. 



45' 



(9.51) 
(9.52) 

(9.53) 

(9.54) 

(9.55) 

(9.56) 
(9.57) 

(9.58) 
(9.59) 



We note from eq. (9.49) 



ds' 



E 

A,B 



AB J 



(9.60) 



Now we define the (anti)-holomorptiic forms 



-AB 



-AB 



^AB,CD^CD 

^(CAB + (^ACecD^DB ± -{^ACeCB - Gac^Cb)) 
1 II . ^ 



= 2^AB =t {UcecB - CBcecA) 



(9.61) 
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At the north pole: 

^12 
^13 
P+ 



'16 



--23 



'24 



-25 



-26 



-34 



-35 



-36 



-45 



-46 



-56 



(ei4 + 623 + «(e24 - eia)) 
(ei5 - 626 + «(e25 + eie)) 

(ei6 + 625 + «(e26 - 615)) 
(623 + ei4 + ^(624 - 613)) 

(624 - ei3 - i(ei4 + 623)) 
(e25 + ei6 + «(e26 - 615)) 
(e26 - ei5 - 'i(e25 + eie)) 

(635 - 646 + «(e45 + 636)) 
(636 + 645 + «(e46 " 635)) 
(645 + 636 + «(e46 " 635)) 



-le 



13' 



-le 



15' 



-le 



13' 



"■=^13' 

4- 

-le 



15' 



-15' 



-le 



35' 



-le 



35' 



-35' 



(9.62) 
(9.63) 

(9.64) 

(9.65) 

(9.66) 

(9.67) 

(9.68) 

(9.69) 

(9.70) 
(9.71) 
(9.72) 

(9.73) 

(9.74) 

(9.75) 
(9.76) 



From eqs. (9.62 1-(9.76 ) we note that only 3 holomorphic forms are independent: 

1 . „ 

6-]^2 ~t~ ^6]^4 I — ^6 



-13 



-15 



-35 



14' 



^15 + *^16 



16' 



^35 + ^^36 



36- 



The line element can also be written as 



where 640 = (6^0)*. At the north pole 



(9.77) 
(9.78) 
(9.79) 

(9.80) 



A<B 



^ (^13^13 + ^15^15 + ^35^35) • 



(9.81) 
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9.2. Review of the construction of CP^ 



Eq. (9.81) reduces exactly to eq. (9.60). 



9.2.2 CP^ as orbit under Spin{5). 

Note that Aab ~ Jabi CL,b = 1 ... 5 generate the Spin(b) subalgebra of Spin{6) 
while = '^'^aG transforms as a vector under Spin{5). 



We define Aq = J7a6 so we can write the projector (9.18) as: 



The projector (9.43) takes the form: 

1 



P^ = -l + lA, + i=(Ai2+A34) 



(9.82) 



(9.83) 



Here we have an extra restriction: under Spin{5) rotations the norm of the vector 



is not affected. From eq. (9.83) we have ^^^a 



4- 



The induced line element is: 



(9.84) 



the constants a, (3 are arbitrary numbers until know. We will come back to this 
point at the end of this chapter. 
Spin{5) rotates Aab as a tensor: 



Cab — RacRbdCcd 



(9.85) 



and it rotates Aq as a vector: 



— RabS,b 



(9.86) 



In analogy to eq.(9.47), d^^^ = —Tr [R ^dR, but now the trace runs over 



the sub-indices a,b = 1 ... 5. For the vector part, dS,a we have: 



-Tr 



R-^dRR-^dRi 



(9.87) 



where represents the matrix of coefficients and ^ ^ the vector ,^J^. Now 
we can use the Maurer-Cartan forms for the group of rotations on 4-dimension: 
R^^dR := —CabLab where L^b are the generators of this representation. For more 
details of these calculations see appendix IfI 
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Finally, the line element at the north pole is: 

as = [ei5 + 625 + 635 + 645] + 2 K^is " ^24) + (cm + 623) J 



(9.88) 



Eq. (|9.49|) reduces to eq. (|9.88|) if we ignore the rotations Cae and we choose 

we recover 



a = 1, P = 1. Then if we choose as a particular case a = /3 in eq. (9. 
the Spin(6) symmetry. 



From eq. (9.88) we obtain geometrical information of CP . We focus into the 



following two points: 

1. Isotropy group of the orbit. 

2. Local form of CP3. 



1) The isotropy group of the orbit 

To find the isotropy group of the orbit, we have to identify the rotations that do 



not affect the projector given by eq. (9.83). Such rotations must be given by the 
orthogonal forms: 



^AB - PAB;CDeCD- 



(9.89) 



'12 



ei2, 



eis = 


1. 

2(^13 


+ 624) 


- P-L 

— 6245 


— 


1. 

2(^14 


- 623) 


- P-L 

— 623, 


^15 — 


1. 

2(^15 


+ 626) 


- P-L 

— 6265 


^16 = 


1. 

2(^16 


- 625) 


- P-L 

— 625, 


p-L - 
634 — 


634, 
1. 

2(635 


+ 646) 


- P-L 

— 645, 


p-L - 
^36 — 


2(636 


- 645) 


- P-L 

— 645, 


^56 = 


656- 







From eqs. (9.90 )-(9.98 ) we have to restrict to rotations on 4-dimensions. The forms 
e;^g, a = 1 • • • 6 are related to rotations exclusive on 5-dimensions. 
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9.2. Review of the construction of CP^ 



From the remaining combination of coefficients we construct the corresponding 
lerators of the isotropy of the orbit, T^, 

ri = i=(A23-Ai4) 

T2 = --^(Ai3 + A24) = 
^(Ai2-A34) 
^(Ai2 + A34) 



= 1,2,3,4 as follows 




-(o-23 - 0-14), 


(9.99) 


^(0"13 + 0-24j, 


(9.100) 


^(C7i2 - 0-34), 


(9.101) 


-((T12 + 0-34). 


(9.102) 



The components do not affect the fiducial projector in eq. (9.83) since: 

[r^,Po]=0, /i = 1,2,3,4. (9.103) 
Ti,i = 1,2,3 fulfil the S'C/(2)-subalgebra 



[Ti,Tj] = lei.kTk. (9.104) 



T4 oc Po generates U{1). 
The orbit turns out to be: 



Spin{5)/[U{1) X SU{2)]. (9.105) 



Note that both prescriptions, developed in sections |9.2.1| and |9.2.2 give the 



same orbit, CP'^, since we can rotate any point of CP^ as SpiniQ) orbit and then 
unrotate it using an element of Spin{5). 

If we follow the construction of CP"^ demanding a S'pin(6)-symmetry we obtain 
a rounded version of CP^. Demanding the less restrictive S'pin(5)-symmetry we 
have a squashed version of CP^. 

2) CP3 is locally of the form x S^. 



As a second observation, in eq. (9.88) we distinguish that the line element is 
composed by two parts: 



A four dimensional part: X]a=i^a5- We identify it as a 5'^-line element. 

t_g,4. 



Then a + (3 is related to the square of the radius of the S^, namely i?^4. It 



can be re-written in terms of the (anti)-holomorphic forms as 

8/3 (e^s^rs + 635635) • 
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A two dimensional part: (ei4 + 623)^ + (eis — 624)^ or in terms of the (anti)- 
holomorphic forms 



4ae+3ei3 



Using the same procedure we used in the calculation of the isotropy group, 



now we look at the combination of indices Cab that appeared in eq. (9i 
We associated with the corresponding combination of generators: 



(ei4 + 623) 

(624 - 613) 



ri = ^/2(Ai4 + A23), 

T2 = \/2(A24 - A13). 



T3 is obtained by demanding that r^, i = 1,2,3 obeys [Tj,rj] = leij^Tk- 
Summarising the above 

ri = ^/2(Ai4 + A23), (9.106) 
T2 = y2(A24-Ai3), (9.107) 
T3 = V2(Ai2+A34), (9.108) 

generate a sphere S*^. 



Then the constant (3 in eq. (9.88) is related to the square of the radius of 
the S^, i?|2- a + P is proportional to the square of the radius of the S^, i?|4 



ds 



n2 
i?2 



R 

[ei5 + 625 + els + els] + ^ [(ei3 - £24)^ + (eu + 623) 



.109) 



In eq. (9.88) the radius of CP^, is 2. Summarising, eq. (9.109) reflects that 
locally CP3 ^S^x S^. 

It is known that CP^ is a fibre bundle with as a fibre and S"^ as a base space 
|51] , in terms of diagrams we have 



^ CP 
i 

The construction of CP'^ as an orbit under Spin{b) allows us to corroborate 
this fact. 
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9.2. Review of the construction of CP'^ 



As we mentioned in the introduction, in order to specify the geometry in fuzzy 
spaces we have to give a prescription for the Laplacian. In chapter [T0| we will get 
the exact relation between radius of the fibre S"^ and a penalisation parameter h 
appearing in the scalar theory of S^. 



Chapter 10 



Decoding the geometry of the 
squashed CP|i 



Ref. |39] presents a prescription for a scalar field theory on a fuzzy 4-spliere. 
However, this theory depends on an additional parameter h and it is defined on a 
Zar^cTQ space, CP3, instead S^. 

The motivation for this chapter is to explore the feature of fuzzy spaces, to 
reflect the geometrical properties, and to explain why the prescription in [3^ works. 



This chapter is divided into two parts: in section 10.1 we present a review of 



the results in [SHj. In section 10.2 we extract the geometrical information from the 



Laplacian via the tensor metric. 



10.1 Scalar field theory on fuzzy 4-sphere 



In this section we will summarise the prescription for working with a scalar field 
on Sp following Ref. |i39j . The action for the scalar field on a rounded CP^ is 



given by eq. ( 10.1 ): 



(10.1) 



where Jab are the Spin{Q) generators defined in chapter 9.2 The constant in 
front of the trace, ^ represents the volume of CP p. 



^CP"^ is a 6-dimensional space. 
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10.1. Scalar field theory on fuzzy 4-sphere 



Since we are interested on retaining just the S'pm(5)-symmetry, we add to 



eq. (10.1) a 5*0(6) non-invariant but SO{5) invariant term given by eq. (10.2): 

= ^Tr^(^[Ja,,<!>]^[Jab<^]-l[JAB<^nJAB<^]^ (10.2) 

SO that we have the overall action for a squashed CP^ 

Sm=Tr{<^Ah<^ + Vm). (10.3) 
The expression for the overall Laplacian, is 

Ah- = ^ I^I[Jab, [Jab, •]] + h{[Jab[Jab, •]] - ^IJab, [Jab, •]])) • (10.4) 
or equivalently 

A, — 

2i?2 

which gives a stable theory for all L if /i G (— 1, oo). 



= ;^fcf(^^ + M2Cf(^^-Cf(^^)) (10.5) 



This form (10.5) is an interpolation between 50(5) and 50(6) Casimirs. As a 
particular case, the Laplacian is proportional to the 50(6) Casimir for h = and 
the 50(5) Casimir for h = 1. 

The probability of any given matrix configuration has the form: 

-Som-hsim 
m = z 



where 



Z = j d[$]e-^[*]-''^^[*] (10.7) 



is the partition function of the model. 

The values of h of interest to us are those large and positive since in the 
quantisation of the theory, following Euclidean functional integral methods, the 
states unrelated to S'^ then become highly improbable. 

Note that we have not specified the potential of the model since the above 
prescription is independent of the potential. The most obvious model to consider 
would be a quartic potential, since this is relevant to the Higgs sector of the 
standard model. 
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10.2 Analysing the geometry encoded in the 
Laplacian 

10.2.1 Mapping the fuzzy Laplacians to the continuum 

The non-commutative product of matrices induces a non-commutative product on 
functions, this is the -k-product. This useful tool allows us to access the continuum 
limit. Let Mi, M2 be two matrices of dimension di, and Mi(^), M2{(,) are the 
corresponding functions obtained by the mapping: 



MiiO := Tr {Vl{OMi) . (10.8) 
Vl{0 is called the projector at L-level and it is contracted taking the L-fold tensor 



product of "P defined in (9.18). Vl{0 carries the coordinates S,. 



The definition of the -k-product is given by 



(Ml k M2) (0 := Tr (^^^(OMiMa j . (10.9) 

For CP^""*^ the k-product can be written as a sequences of derivatives on the 
coordinates ^. For our proposes we will just follow the prescription given in Ref. 



18] . Now we present the formal definitions of the set of coordinates ^ in eq. (9.18): 



Ub := Tr {Vl^ab) , 



(10.10) 



where A^^ are proportional to J'ab as in eq. (9.14). Note eq. (10.10) is consistent 



with eq. (9.18) via eq. (9.16). 



A first interesting fact is that the commutator of J^ab maps into the covariant 
derivative: 



CabM{0 ■■= Tr{^ 



Jab,M 



V2 



JAB;CDdcDM{i). 



The quadratic Casimir operators are defined as: 



[Jab, [Jab-, •]] 

[Jab, [Jab, •]] 



a 
a 



50{6) 
2 

50{5) 



(10.11) 
(10.12) 



(10.13) 
(10.14) 
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10.2. Analysing the geometry encoded in the Laplacian 



What we want to calculate is the image of A/i, eq. (10.5 ), under the ^-product 



map. First we define the image of eqs. ( 10.13 )-( 10.13 ) in eqs. ( 10.15 )-( 10.16 ): 



Jab, [Jab, M] 

Jab,[Jab,M] 



we can rewrite as 



where 



Kg 



K5 



JAB,CDdcD {JAB,EFdEF) 

PcD;EFdcDdEF + JAB;CD{dcDJAB;EF)dEF 
J ah, C ode D {Jab,EFdEF) 



^10.15) 
^10.16) 

^10.17) 
^10.18) 



(10.19) 
(10.20) 
(10.21) 



A very useful expression to simplify the calculations is the contraction of the 
complex structure to the partial derivative 



Jab-cdQcd = V2{iAcdcB - iBcdcA) ■ 



(10.22) 



We are interested in extracting the metric tensor Q comparing the relevant con- 
tinuous Laplacian with the general form: 



(10.23) 



G^'^d^d, + {d^Cn du + ^G'^' (s^Vg) d,. (10.24) 

V G ^ ^ 



For the case when we retain the full 5'pm(6)-symmetry we have the Laplacian 
C^^\ The associated metric tensor is just Pab;CD as it can be verify from a 
straightforward calculation for kq: 

1 



= Ij^^AB - '^iAciBodABdcD - 8S,ABdAB- 



(10.25) 

For the 5'pin(5)-symmetry case the corresponding image to the Laplacian 



eq. (10.4) is defined proportional to £| 
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Tr ( AhM 



2E? 



'-h 



then 



C\ = KQ + h (2k5 - kq) . 



(10.26) 



(10.27) 



Our guess for the tensor metric related to C\ can be decomposed into a pure 
S'pin(6)-symmetry part (given invariant part denoted 

by as in eq. (|l0.28|): 



qAB-CD _ pAB:CD _|_ i^j^AB;CD 

X is related to (2k5 — Kg), the term that breaks the Spin{6) symmetry. 
After a straightforward calculation we get: 

1 



(10.28) 



2k5 - K6 = ^dlf, - '^iabicddacdbd " ^Ueibedacdbc - ^abdab 



(10.29) 

The tensor x^^'*-"^ is obtained comparing the terms with second derivatives in 



eq. ( |10.29| ) to G^^'^^dAsdcD 

j^ab;cd ^^ab;cd 2 /' ^"^c^fed ^ad^i 



be 



a6:cd 



2 ^(J'^c^eb^ed _ ^ad^6b^6c _|_ ^bd^Ga^Qc _ jfec^6a^6d^ ^-^q ^q-j 
= X"^''^^ = X"^''^^ (10.31) 



We can rewrite X as: 



j^ab;cd pab;cd _ j\jab;cd 

2 



ji^ab;cd ^ j jac^6f)^6c( ^ad^db^dc _j_ ^6d^6a^6c ^bc^6a^6d 



with 



Some properties of M"^'^'^ and X°'^''^'^ are: 



-ji^ab;cd jy^cd^ab jy^ba-jCd 



X 



ab:cd 



j^cd;ab j^ba;cd 



X^ = X. 



(10.32) 



(10.33) 



(10.34) 
(10.35) 
(10.36) 
(10.37) 
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Now the traces amount to 



j^ab;ab 
j^ab;ab 



1 



5 + 2^"^^'"^ - lec^ 



6bc6b 



2. 



(10.38) 
(10.39) 



Notice X"'^''^'^ is a rank-2 projector while M°''''^'^is a rank-4 projector. 

In order to invert the tensor metric we calculate the products between projec- 
tors pAB'CD^ ^afe;cd ^ab;cd. 



pab;ef pefyCd 
pab;EF -^EF,cd 
pa6;EF j^EF,cd 
pab;EF j^EF,cd 
pa6;EF j^EF,cd 



pab;ef _ _j^ef,cd 

4 ' 

pab;ef j^efyCd j^ab;cd 



pa6;ef j^ef, 



cd 



0, 

1 



pab]ef j^ef,cd _j]^ab]cd 

~ 2 

pa6',ef j^ef^cd pa6;cd 

0. 



The metric fulfils: 

Now we assume the covariant tensor metric to take the form 

GaB;CD = PaB;CD + OiXAB,CD, 

where a is a constant to determine. 



From eq. (10.46) we distinguish two cases: 

Gab\CD ^ 



= {P + hXr-^^^{P + aX)cD;ef 

= P'^^'^f + {ah + a + h)X''^'^f 

Comparing eq. (10.48) to eq. (10.46) we gelj^a = 

„ „ h 



'ab;cd 



ab:cd 



1 + 



X 



ab,cd' 



(10.40) 
(10.41) 
(10.42) 
(10.43) 
(10.44) 
(10.45) 

(10.46) 
(10.47) 



(10.48) 



(10.49) 



^We have to demand ft, 7^ —1 in order to invert the tensor covariant metric tensor 

qAB-CD 
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pa6;ef 



(10.50) 



From eq. (10.50) we observe that for this case GcD;ef does not depend on h. 



=^ Ga6;cd — Pa6;cd 



(10.51) 



Finally we summarise the results in eqs. ( 10.49 )-( 10.51 ) in eq. (10.52): 



Gab ,00 = Pab:CD 



h + l 



Xab,cd- 



(10.52) 



Observations 



Eq. (10.45) suggested that Xab-cd and Mab-cd project onto separated spaces. As we 
demonstrate in section 9.2.2 CP^ is locally of the form S"^ x S"^. 



Xab;cd projects onto a 2-dimensional space that should be the fibre 5^ while 
Mab-cd projects onto a 4-dimensional that should be the base S^. This information 
can be precisely obtained analysing the contributions to line element ds^ due to 
the tensors Xab-cd and Mab-cd- This will be achieved in the subsequent section 
110.2.21 

Let us define fl, the symplectic structure as: 



CD 



^AB,EF = GaB,CdJeF 



(10.53) 



It is clear when h = (full 5]5in(6)-symmetry case) ^ab,cd = Jab;CD- In general 
we have 



^ab,c6 — Pab,c6, (10.54) 
^ab-,cd = {l + h)Pab-,cd- (10.55) 

10.2.2 The induced line element ds^. 

The line element, ds^ at the north pole is defined as follows: 

ds' := G%.^cDe''''e^^, (10.56) 
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10.2. Analysing the geometry encoded in the Laplacian 



where e^^ are defined as the coefficients of infinitesimal rotations under Spin{6), 
i.e. the Maurer-Cartan forms, R^^dR = —cabLab-, where Lab are generators for 
rotations in 5 dimensions. 



Note that in the case /i = (full 5'pin(6)-symmetry) eq. (10.56) reduces to 

(10.57) 



ds' = PAB;CDe^''e^^ 



as it was expected. 

At the north pole we have 



c5^2 



2(6^=^) 



vO pib cd 



,13 „24\2 



(10.58) 
(10.59) 



From eqs. ( 10.58 )-( 10.59) we verify that Mab;cd projects to the base space 5^ 
and Mab-cd projects to the fibre S"^. The line element at the north pole depending 
on the parameter h is 



1 



Worthpole = (e?5 + 625 + 6^5 + 645) + (ei3 " 624) + (cu + 623) 

(10.60) 

We can compared eq. (10.60) to eq. (9.88). We identify the radius of S"^ and 5^ as 



F?2 

i?2 



1 + 
1. 



^10.61) 
^10.62) 



As we saw in section 10.1 h takes values in the interval (— l,oo). We analyse 



some particular cases of eq. (10.60): 



1. /i = 0. We recover the Spin{Q) symmetry. The radius of the rounded CP'^ 
is i?2 = 2. 



2. h 



oo. In this limit 



3i 

i?2 







then, the limit /i — > oo corresponds to shrinking the fibres to zero size. 
The radius of remains finite: 



R^ 



1. 
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3. h 



-1. Note that h cannot take the exact value 1 since for h = 1 the 



tensor metric 

qAB;CD 

in eq. (^10.28^ cannot be inverted. For this hmit we 

have 



00, 



P2 



fl2 - ' R2 

This Hmit corresponds to making the size of the fibre infinitely large. 



Chapter 11 



Conclusions from part III 



We reviewed the prescription given in Ref. [39] of the scalar field theory on 
a fuzzy 4-sphere. 

Since 5^ is not a phase space the construction of its fuzzy version has some 
complications. The construction starts defining the matrix coordinates Xa 
as 

'''' - vw+Ar' 

where Ja belong to the (^j ^) representation of Spin(5). They are matrices 
of dimension (11 = with L an integer number. From JaJa = 

L{L + 4)1 follows an equation which holds for the matrix coordinates of a 
4-sphere 

5 

XaXa = R^l. 

a=l 

In the limit L — > oo the matrix coordinates commute, then we have a 
matrix approximation to at algebraic level. 

The complications emerge for a finite L where the five matrix coordinates 
Xa, a = 1, . . . , 5 do not provide a complete basis for the algebra of functions. 
To circumvent this problem we have two options: 

1. To include more generators to complete the basis. This leads us to 
include more coefficients in the expansion of a function that are not 
related to degrees of freedom of S"^, i.e. extra degrees of freedom. Then 
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the constructed space is not Sp, it turns out to be CP^. To recover 
a scalar theory on Sf,. We implemented a penalisation method for all 
the non-Sp modes in CP^. 

2. A second alternative consists of projecting out the non-S'^ modes. This 
leads us to deal with a non- associative algebra. We found that following 
the previous alternative we got the projector of the non-S^ modes 
required in this second option. 

The penalisation method introduced in the previous point consists of the 
following: 

— We start defining an initial action on CP|., iSof^]. 

Then we modify So[^] by adding a term S'/[$] gives a positive 

value for those field configurations associated to the non-S'^ modes and 
it is zero for those of S'^. 

— Then we construct the overall action as 

= 5o[$] + ^S7[$] (11.1) 

where /i is a penalisation parameter in the interval (—1, oo). 

— The probability of any given matrix configuration has the form: 

-Sol^]-hSil^] 

n^] = — z — ^^^-^^ 

where Z = J d[^]e~^^^^~^^^^^^ is the partition function. For h — > oo 
the states unrelated to become highly improbable. 

The resulting action is 

S[^] = Tr{^Ah^ + Vm), (11.3) 

where 

A. = ^(cf(^) + M2Cf(^^-Cf(^^)) (11.4) 



is written in terms of quadratic Casimir operator of the groups 5*0(5) and 
SO{t 
CP^ 



SO{6), C2'^*'^'' and C2*^'-^^ respectively. is the square of the radius of 
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Since the defined action S[^] for h large and positive contains only modes, 
the resulting model describes a scalar field theory for 



F- 



The fuzzy spaces can well retain the geometrical properties of the discretised 
space. This nice feature allowed us to provide a exact geometrical interpreta- 
tion why the penalisation procedure described above works. The Laplacian 



is given in eq. (11.4). 



CP'^ can be constructed either as a Spin{6) orbit or as a Spin{5) orbit. The 



Laplacian in eq. (11.4) corresponds to CPp with Spin{5) symmetry. We 
restore the Spin{6) symmetry in a particular case of eq. ( |11.4 ) with h = 0. 



It is known in the literature that CP'^ is a non-trivial fibre bundle over 
with as the fibre. The construction of CP|, as Spin{5) orbit clarified this 
fact. 

Using coherent states techniques we extracted the covariant tensor metric 



in eq. (11.4). Once we inverted the tensor metric we found the line element 
ds'^ in terms of the parameter h as 

5 5 

ds^ = YTh ^ Xab;cde'''e'^'' + - Yl Mab;cde'''e''', (11.5) 

a,b,c,d=l a,b,c,d=l 

where e""'^ are the Maurer-Cartan forms for Spin{5). Xab;cd is a rank-2 
projector and Mab-cd is a rank-4 projector. 

Xab-cd projects to the fibre S"^ and Mab-cd projects to the base space S^. The 
line element in eq. (|ll.5l) shows that locally CP^ = S^xS'^. 



We identified the radius of the fibre and base space as 

-R^2 1 



i?2 
i?2 



1 + 
1. 



(11.6) 
(11.7) 



As h 



oo we have 



0. 



The meaning is that the S"^ fibres shrink to zero size while the radius of 
remains finite. 
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The suppression of the non-S*^ states in the field theory on CP^ corresponds 
to a Kaluza-Klein type reduction of CP'^ to S^. 



As we mentioned in section 9.1 following our construction of CP^ as a 
Spin{5) orbit we are able to give the prescription for the projector to 5^ 
modes 



-ps' = nni 



n p,SO(5) , 
^2 ~ '^r. 



_p50(6) _ . 
n=lm=l 2^2 ^n,m 

L n „ 

1 + n n c^o,.,"' , ■ 

n=l m=l ^2 ■^'^n,m 



Chapter 12 

General conclusions and 
perspectives 

In this thesis we have presented a study of scalar field theory on fuzzy spaces. 
The main motivation of our work was to explore the fuzzy approach as a non- 
perturbative regularisation method of Quantum Field Theories. For this pur- 
pose we chose a hermitian scalar field theory in a three dimensional space, with 
Xcj)"^ potential. This model is perturbatively super-renormalisable (i.e. there are 
only a finite number of perturbative diagrams that require renormalisation) . Our 
non-perturbative regularisation consisted of a fuzzy two sphere for space and 
a lattice for Euclidean time. We performed Monte Carlo simulations and ob- 
tained the phase diagram of the model. The study of such model via a stan- 
dard lattice regularisation leads to a phase diagram consists of a disordered and 
a uniformly ordered phase separated by a continuous second order phase tran- 
sition that is governed by the Ising universality class. In contrast to the stan- 
dard lattice regularisation [91, in this new model we found three phases, two 
are the disordered and uniformly ordered phases but they are separated by a 
third new phase of non-uniform ordering. This third phase is a property of the 
non-commutativity of the regularised model and has arisen in other studies in 
the literature and has variously been called a striped phase (Gubser and Shondi 
|26] . Ambj0rn and Catterall [28], Bietenholz et al. [35]-[3^) or a matrix phase by 
Martin [30]. We find that the three phases meet at a triple point characterised 
by (Ar,mV) = ((41.91 ± 15) Af-0-64±0-20ij-i.28±o.25^ _(12.7 ± i)i?-i-92±o.20)^ 
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section [7] which inevitably runs to the origin as the matrix size is sent to infinity. 
The imphcation is that the model we have studied in the end does not capture 
the non-perturbative physics of the flat-space Xcp^ field theory. This was not un- 
expected since perturbative studies [29] found that this fuzzy sphere model suffers 
from Ultraviolet-Infrared (UV/IR) mixing, and the non-uniformly ordered phase 
is a non-perturbative manifestation of this phenomenon in the neighbourhood of 
the ordering transition. It is conjectured (but still has not been demonstrated 
in numerical studies) that the introduction of another term into the action will 
suppress UV/IR mixing and bring the model into the Ising Universality class. We 
have not pursued this issue in this thesis as the introduction of this additional 
term involves an additional parameter in the phase diagram and it was necessary 
to proceed in steps. First it was essential to understand the phase structure of this 
three parameter model before studying how the phase diagram is deformed by the 
additional parameters. 

The new model is naturally a non-commutative model and is of interest also 
as it is a non-perturbative regularization of a non-commutative field theory, that 
of a scalar field theory on the Moyal-Groenewold plane. The new phase seems to 
be a characteristic feature of such non- commutative theories. On the fuzzy sphere 
it is characterise by the dominance of non zero angular momenta (i.e. / > 0) in 
the ground state of the model and implies that even though rotational invariance 
is preserved in the regularisation process it is spontaneously broken by the ground 
state. 

An advantage of the study is also that different limiting spaces (e.g. the com- 
mutative sphere, the Moyal plane and the commutative flat space) can be obtained 
scaling i? as a power on N and taking N — > oo. Our analysis shows that in all 
the limits we have considered so far the phase of non- uniform ordering survives. 
In particular (as mentioned above) the commutative flat model is not in the Ising 
universality class. However, we expect that a normal ordering in the vertex, i.e. the 
introduction of a counter-term to cancel the tadpole diagrams, would return the 
model to the Ising universality class. This modification in the action is equivalent 
to reinforcing the fuzzy kinetic term and can be achieved more simply adding a 
term <^>{C'^)'^^/{A'^R'^) where A is a momentum cutoff to the action. An analogous 
prescription of adding higher derivative terms to the quadratic terms should be 
applied to all fuzzy models, in such a manner that all diagrams are rendered finite. 
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if the commutative theory is required in the infinite matrix hmit. 

This three dimensional model is amenable to an alternative treatment as a 
Hamiltonian on the fuzzy sphere where the Hamiltonian would be 

47ri?2_ /1„, , 1 , . 1 2^2 . A 



n = Tr -n^ + — ^^>£^$ + -m^$^ + -^>' 

\2 2R^ 2 4 

with n = 6. 

At strong coupling A it is expected that the model effectively depends on 
less parameters. We found that our simulations in that regime reproduce the 
predictions for the chain of matrices |17] where the contribution of the fuzzy kinetic 
term is neglected. This model is interesting from the String Theory point of view, 
it has been studied in Refs. [6lj- [65] where it is known as c = 1 model. 

At the technical level, the simulations of scalar theories on fuzzy spaces are 
slower than in their lattice counterparts since the fuzzy models are intrinsically 
non-local. Nevertheless we have found that the simulations quickly converge as 
the matrix size is increased and we can capture the characteristic behaviour of the 
observables at very small matrix size. It is, however, expected that true advantages 
of the fuzzy approach only emerges with the simulations involve a fermionic sector 
or the models of interest are supersymmetric. 

A second, theoretical aspect of the thesis was the presentation of four dimen- 
sional models on a round fuzzy approximation to 5^, via a Kaluza-Klein reduction 
of CP^. Though we did not perform numerical studies in this case, from the stud- 
ies that we have done so far, we are in a position to conjecture the structure of 
the phase diagram for the CP^ and S*^ models. Since, the disorder ed-non- uniform 
ordered transition line seems to be universal in the class of models where space- 
time is modelled by a fuzzy space, we expect it to arise in theses models also. It 
corresponds to the dominance of the potential term and the models become pure 
potential one matrix models (with, in our case, a cj)'^ potential). The large A^ limit 
of these models are solved in terms of the density of eigenvalues. This density 
undergoes a transition from a connected to a disconnected density as the potential 
well is deepened by making the mass parameter, m^, more negative. It is this 
separation of the eigenvalue spectrum that occurs at the disorder ed-non- uniform 



order transition line. In the pure potential model (see eq. (6.50)) the transition 
curve is given by be = —2\/Nc where h is the total coefficient of the quadratic po- 
tential term and c that of the quartic term. This disordered phase should give rise 
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to a uniform ordered phase as the mass parameter is furthered decreased. Again 
all these models should exhibit UV/IR mixing. The study of adds further com- 
plications due to the need to squash CP^. A further complication is the need to 
introduce the additional UV/IR suppressing term. However the principal conclu- 
sion of this thesis is that we see no insurmountable difficulty to the implementation 
of the scheme implemented in this thesis and the extensions outlined above as a 
regularisation scheme for quantum field theories. It has the distinct advantage 
that it is also a natural regular method for non-commutative field theories. This 
constitutes an interesting numerical experiment for the Kaluza-Klein reductions 
in the fuzzy context. 

Although there is still a lot to do to reach the high acceptance status of Lattices 
Field Theories, to the date the fuzzy approach have overcome the very first tests. 
We may take the present work as a starting point to continue the exploration of 
this possible alternative in future studies on fuzzy spaces. 



Appendix A 

A small description of the 
Monte Carlo method 



Let ^ he a configuration of tlie relevant field. The probability for to this configu- 
ration is given by 

e-5[l>] 

m = (A.i) 

where S[^] is the Euclidean action of the system in the configuration ^ is 
called the partition function, 

Z = J D[<^]e-^^'^\ (A.2) 

where / -D[$] denotes the integration over all field configurations. The expectation 
value of the observable O is define by the expression: 

(0) = J nm^Om. (A.3) 

The idea of the Monte Carlo method is to produce a sequence of configurations 
{^i}, i = 1, 2, T/v/cj^ and evaluate the average of the observables over that set 
of configurations. In this way the expectation value is approximated as 

{0)^ — J2 Oi, (A.4) 
Tmc ^ 

where Oi is the value of the observable O evaluated in the i-sampled configuration, 
i.e. Oi = 0[$i]. 

^Tmc- Monte Carlo time 
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A.l. The Metropolis algorithm 



The sequence of configurations obtained by Monte Carlo have to be represen- 
tative of the configuration space at the given parameters. 

A.l The Metropolis algorithm 

The concrete way that the Metropolis algorithm works is the following: 

• Start with a configuration ^mit = and generate another configuration 

'^test- 

• Compute AS := S[^test] - S[<^o]. 

— If AS < then ^test is accepted, i.e: 

^1 = '^test- 

— Otherwise the number e~^'^ is compared to a random number, ran S 
[0, 1]; if e~^^ > ran the configuration ^test is accepted, in the other 
case it is rejected, i.e. = ^q. 

• Set ^init = ^1 and compare it with another configuration ^test- 

Before we measure it is necessary to perform a number of steps, in order to obtain 
stable values for the observables of interest. This is called thermalisation time. In 
our study the thermalisation time was estimated from the history of the actioiij^ 

The standard way to propose the following configuration 

The variation of the configuration <I>(t) is performed element by element, 

^{t)ij ^'{t)ij = + aij (A.5) 

where aij G C is a random numberj^ Its real and imaginary part are in the interval 

and the chosen so that we preserve hermiticity of 



the field, i.e. 



^ The picture obtained plotting the Monte Carlo step i vs. the action at the configu- 
ration S[<^i\. 

^In general the real and imaginary part are different random numbers. 
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One Monte Carlo step correspond to updating all entries of each t 
1, . . . , Nt sequentially once. 



A. 1.1 Modifying the Metropolis algorithm 

The standard Metropolis algorithm works well when the observable fluctuates 
around one value. In our study we found that the value of the observable fluctuates 
around several values]^ This was interpreted in the sense that the effective poten- 
tial of the system has several local minima. It was also observed that for certain 
values of the parameters there is no tunnelling between the different minima for a 
very long history. Then the result depends on which minimum the system falls in, 
and this typically strongly depends on the starting conditions. 



Independent simulations 

As a first try to handle this dependence on the starting conditions we perform 



many independent simulations. The algorithm described in A.l changes as follow: 

• Divide the T/v/c-Metropolis steps into sim parts. 

• Choose a starting configuration and thermalise. 



Perform the loop describe in A.l TMc/sim times. 



• Repeat the previous steps sim times until you collect Tj\,fc'-configurations. 

This method is useful to check if the results do not depend on the starting 
conditions, but in some cases this method fails because the expectation values of 
the observables depend on the way that the sim different initial conditions are 
chosen. 

Let us characterise each minimum by the expectation value of the energy. 
In figure 5.17 of chapter |5] we showed that there is a clear difference between 



results with cold start vs. hot start. Between simulations with a hot start no 



tunneling is observed, as figure A.l illustrates. 

We call it multilevel behaviour, some examples of it appears in appendix jcj 
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A.l. The Metropolis algorithm 
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Figure A.l: History of the action for 6 different hot starting conditions at 
m2 = -24, A = 1.25, ^ = 8, = 12. We take a new start every 500, 000 
Monte Carlo steps. 



In figure A.l we observe the history of the energy for 6 independent simulations. 



We observe 4 different mean values in the sense of (A.4). This indicates that the 
space of configurations is divided into separated subspaces since we do not observe 
tunnelling in the same simulation. 

According to the mean value of the energy, in the second and third simulation 
the system falls in the same minimum, characterised by the central value {S) ~ 
—520. We have a similar situation for the fourth and sixth simulation where the 
expectation value of the energy is (S) w —622.5. 

If a history begins with a random hot start, the probability to be trapped for a 
very long period in a specific action minumum is what we denoted as the "size of 
the minimum" in Chapter [5} In this situation it is obviously problematic to rely 
on a single run, or on the average over a few runs. 



A. 1.2 An adaptive method for independent simula- 
tions 

A second attempt to sample the configurations takes into account the size of the 
minima. To do this we stick the independent simulations in the following way: 

• Divide the T/vfc-Metropolis steps into sim parts 
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Choose a starting configuration and then do the loop described in section 
A.l for TMc/sim times. Keep the last configuration, this will be the ^mit 



for the next loop. 



• Choose a new starting configuration and thermalize to get ^test, and test it 
as in 



A.l with ^init obtained from the previous independent simulation. 



Figure A. 2 shows the Specific Heat for N = 12, X = 22/12, R = 4 obtained 
by the three different methods of measurement discussed in this appendix. 
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Figure A. 2: The Specific Heat measured by the different methods at A = 1.83, 
i? = 4, A = 12. 



As we can observe, different methods give different results and this is because the 
value of the observables strongly depends on the way of measuring. As an example 



of the statement above we present in figure A. 3 the histogram for the action at 
m? = —1.7 of the figure A. 2 
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A.l. The Metropolis algorithm 




Figure A. 4 shows the histogram for Tr(<&) with the same setting as in figure 



12 



10 - 



6 - 



4 - 



2 - 







"1 1 1 1 H 

Independent simulations - 
Adaptive Metropolis - 



± 



]\ L 



_L 



_L 



-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 

rr($) 

Figure A.4: Histograms for Tr($) at = -1.7, A = 1.83, ^ = 4, = 12. 



The thermalisation problems reflect the fact that there exist a large potential 
barrier between the different subspaces. As consequence we cannot sweep all the 
space of configurations, therefore the expectation values of the observables are not 
reliable. 
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A. 2 Methods to estimate the error 

In this section we give a brief explanation of the methods we used to estimate the 
statistical errors in om' simulations. 

We start by defining the mean O over a subset of {Oi}J^^^ of n elements 
(n < Tmc) 



1 

n ^-^ 



i=l 



(A.6) 



The standard deviation is defined as 



n — 1 



n — 1 



02 _ (0)2 



(A.7) 



If the samples are statistically independent eq. (A.7) gives a good estimation of 



the error. If this is not fulfilled, then the correct expression for a is (see Ref. |52]): 



1 + 



a 



2t 
An 



n — 1 



02 - (0)2 



(A.^ 



where r is the autocorrelation time and An is the Monte Carlo time interval at 
which the samples were taken. It is related to the total number of samples by the 
expression n = 



For large n and 2r 3> An, eq. (A. 8) turns into eq. (A.9|: 



2t 

Tmc 



02 -(0)2). 



(A.9) 



Now the problem is to estimate the autocorrelation time. The formal expres- 
sion for the autocorrelation time (see Ref. [53]) is 



k=l ^ ^ 



(A.IO) 



with 



A{k) = ^ 



(0,0,+fc)- (0,)(0,-) 



(A.ll) 
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A. 2. Methods to estimate the error 



A. 2.1 Binning method 

This method is also called blocking method. The idea behind this method is 
to divide the vector of measurements {Oi}J^'^ into rib blocks (also called bin 
number), then we evaluate the observable for each block to obtain a new vector 
of measurements {Oj}^^^. The error is estimated as if these new measurements 
were statistically independent, then it is obtained via eq. (A.7| replacing n by n^. 
The principal disadvantage of this method is that the error may strongly depends 
on the choice of rib. Then one should test several values of rib and keep the one 
where the error is maximal. 



A. 2. 2 Jackknife method 



This procedure can be consider a re-sampling method. It also starts by dividing the 
vector of measurements into rib blocks. Then one forms Nb large blocks cointaining 
all data but one of the previous binning blocks 



O 



jackknife 



Tmc -O-k-Oj 

Tmc - 3 



(A.12) 



where Oj is the average in the j block, k is the number of samples in each block. 
The error is then calculated as follows: 



a 



Nb-1 



\ Nb 



Yjackknife 



o 



(A.13) 



A. 2. 3 Sokal-Madras method 

This method is based on the estimation of the autocorrelation time given by 



eq. (A. 10), see |54|. The error is given by eq. (A. 9). In the case that the au- 
tocorrelation time turns out to be below 0.5 (the samples are decorrelated) , we 



take the standard error given by eq. (A. 7). 



In our simulations we compare the error given by these three methods. In the 



case of the observables (/?^,,, ipQ^ xi and E - eqs. (4.27), (4.30), (4.32) and (4.34) 



respectively ~ the errors given by the three methods are compatible. 

For the Specific Heat (or the susceptibilities of the different modes) the most 
careful estimate was given by the Sokal-Madras method (largest errors). 
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A. 3 Technical notes 

The runs were performed on two different clusters: 

• Berlin cluster. 

- 450MHz Athlon (pha4-pha9) 

- 900MHz (phal5-phal9) 

- 2.66GHz P-4 (pha30-pha33). 

• Dublin cluster 

- 3.06GHz Intel Xeons-dual (cluster0-clusterl5) 

- 2.80GHz Intel Xeon-dual (Gibbs) 

- l.SGhz P-4 (schrodinger, hamilton, oraifertaigh, lanczos) 

The code is written in C + + and uses the Message Passing Interface (MPI) . 



Appendix B 

Polarisation tensors for SU{2) 



In this section we want to describe some generalities of the polarisation tensors 



Yi^m, which form a basis for the matrices $ in eq. (4.13 ). They are the eigenvectors 



of the operator defined in eq. (3.9) -the fuzzy version of the angular momentum 
operator £^ 



i=l 

where Lj S Mati+i- So we have the set {lim}«<L,m<« such that 



C^Yi^ = l{l + l)Yi. 



Imi 



(B.l) 



(B.2) 



Following Ref . [55] , their algebra is 



^imi Yl2m2 



L + 1 



5](-l)^+V(2^i + l)(2/2 + l) 
f h h I' 



I'm 



L/2 L/2 L/2 



~,l'm' 
-'hmil2m2 



(B.3) 



Yvm' (B.4) 



> are the Winner 3mi- symbols -see Ref. 1561- I57i - and C ™ 
L/2 L/2 L/2 J 

are the Clebsch-Gordan coefficients. 



hm\l2m2 
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B.l. Explicit form of the generators of SU{2) IRR of dimension (11 



Their normahsation chosen as: 



An 
L + 1 



T^iYLyi'm') =Su'5mm'- (B.5) 



B.l Explicit form of the generators of SU(2) 
IRR of dimension di 



As we mentioned in section 3.2 Lj are the generators of the SU{2) irreducible 



representations (IRR) of dimension N = L+1. They satisfy the Casimir constraint: 

3 1 

Clui2) := E^^^' = - 1) • 1 • (B.6) 

1=1 

There is a weU established procedure to construct these generators in an 
arbitrary representation. It operates in the Cartan basis where we work with 
:= Li + iL2, := Li — iL2 and := L3. 

[LJ« = 1'"''^'-^*' (B.7) 

[ otherwise 

I otherwise 

[LA., = [ ^'^''-'^ (B.9) 
otherwise 



Note that (L+)t = L_. 

The commutation relations read: 



[L+,L_] = 2L:, . (B.IO) 



^In fact the operators defined in [SS], T;™, are essentially our polarisation tensors up 
to a factor: 



B Polarisation tensors for SU{2) 
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B.2 Explicit construction of the polarisation 
tensors 

l/m can be constructed as traceless polynomials of order / on Lj, for example, for 
/ = 0, loo oc 1. Demanding j^Tr {YqoYqq^ = 1, we arrive at loo = (also 
see Ref. [58]). 

For / = 1 we have Yio oc L^, Yi^i oc L-|_ and li,-i oc L_. The normalisa- 
tion can be found calculating the trace of the squared of L^, L^^-, L^. Following 



our definitions (|B.7l)- (|B.9|) we found Tr ( J = M^i±li(^±?)^ Tr ( L^, L+ 



Tr (L„L+) = _ Then we choose Yi +i = e*"?^ J^^J=L+ and so on 



for the remaining m's. The phase e*"^ has to be fixed demanding eq. (B.3) to hold. 
Summarising: 

Yoo = -^1; (B.ll) 
V47r 



Yi+i = tj- , ^ L+ (B.12) 

fio = \- , ^ =i^z (B.13) 

Yi_i = iW-?^ , ^ (B.14) 

V27r/L(LT2) ^ ^ 
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Aside results 



C.l Criteria to determine the phase transi- 
tion 



To sketch the phase diagram we compare the fohowing two criteria: 

• The criteria of the susceptibihties. The first phase diagram was revealed by 
searching for the values of where the two-point function of the zero mode 
has its peak for a given value of A. 

• Specific Heat Criterion. A second phase diagram was found searching for 
the values of where the specific heat peaks. 

We found that for values of A < for R fixed both criteria roughly coincide, as 

16, A = 0.44, R = 4: 



it can be seen in the figure |C.l| for N 
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Figure C.l: Susceptibilities Xo and Xi and the specific heat at A 
R = AN = lQ. 
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C.l. Criteria to determine the phase transition 



We observe that the specific heat an the zero mode susceptibihty xo in eq. (4.31 ) 
roughly peak at the same value of m^. This can be explainrd analysing the partial 
contributions to the action: 
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Figure C.2: Partial contributions given by eqs. (4.40)-(4.43) to the internal 



energy eq. (4.34) at A = 0.44, R = A, N = 16 



We observe that the main contribution comes from the kinetics terms. The 
kinetic fuzzy term selects the configuration where the zero mode is leading. 
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Figure C.3: Comparison between the critical points obtained by the critera of 



the Specific Heat -eq. (4.35)- and of the zero mode susceptibility in eq. (4.31 ), 
at -R = 8, = 12. We observe that the critical points overlap within the 
error bars. 
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C.2 Free field results 



In the case A = 0, following Ref. the expression for the space correlator (4.22) 
for / = is: 



(cSo(O)coo(O)) := {ifl) 



1 



(C.l) 



Figure C.4 shows that our simulation results agree with this formula. 
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Figure C.4: {iff) for A = 0, = 4, = 12. 



Appendix D 



Tables 



Table D.l is organised as follows: 



1. For each pair (N, R) it gives the maximal value of A simulated — for several 
values of m? — free of thermalisation problems around the phase transition. 

2. Next it indicates the critical value, m^. This is done with the purpose of 
having a reference of the magnitude of the simulated values of . 

3. Then it contains the corresponding parameters A,B,C and D defined in 



eqs. ( [7l3l )-( [7l5| ) (see chapter [7] for more details). 
4. Finally it present the proposed fit for the transition curve I-III of the form 



(6.45) (in some cases it indicates the range of A for the proposed fit). 
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Tabic D.l: Critical values for the disorder ed-ordered non- uniform phase tran- 
sition 



Maximal critical 
parameters 


Proposed fit for 
the transition curve 


TV = 8 R = A 7VA= 1300 m2 = -24±0.5 
A = 2tt D= 100.531 Be = -2412.74 C = 65345.1 


TV?e = -0.099(ArA)°-^6 


N = 8 R = 8 NX=2500 m2 = -13.8±l 
A = 2n D = 402.124 = -5549.31 C = 502655 


= -0.065(iVA)0-69 


N = 8 R=16 NX = 5000 m = -9 ± 0.3 
A = 2n D = 1608.5 Be = -14476.5 C = 4.02124 x 10^ 


m-^e = -0.026(7VA)°-69 


N=12 i? = 4 iVA= 1200 rn2 = -24 ±0.5 
A = 2tt D = 100.531 Be = -2412.74 C 60318.6 


m^e = -0.15(A/-A)0-^i 
for NX > 400 


TV =12 R = 8 7VA = 2000 = -16.5 ± 0.5 
A = 2n D = 402.124 B^ = -6635.04 C = 402124 


m^e = -0.05(A^A)°-^6 
for A^A G [600,2000] 


N=12 R=16 NX= 1800 = -6 ± 0.2 
A = 2'K D = 1608.5 Be = -9650.97 C = 1.44765 x 10^ 


= -0.037(AfA)°-68 
for A/^A > 300 


AT = 12 R = 32 NX= 1600 = -2.2 ± 0.2 
A = 27r D = 6433.98 Be = -14154.8 C = 5.14719 x 10^ 


m'^e = -0.016(7VA)0-6<5 

NX > 20 


iV=i2 i? = 64 NX = 60 = -0.1125 ±0.006 
A = 2tt D = 25735.9 Be = -2895.29 C = 772078 


= -0.013(A/'A)0-5i 
for A^A > 15 


= 12 i? = 100 NX = 30 = -0.0535 ± 0.0025 
A = 2n D = 62831.9 Be = -3361.5 C = 942478 


to2^ = -0.005(AfA)°-64 
for 7VA > 5 


N =16 R = 2 NX=120 fh? = -6.8 ± 0.25 
A = 2-K D = 25.1327 B = -170.903 C = 1507.96 


m'^e = -0.155(AfA)0-^9 


N =16 R=A NX= 1200 m=-27±0.5 
A = 27r D= 100.531 B = -2714.34 C = 60318.6 


m2^ = -0.084(ArA)0-Si 
for A/'A > 400 


A^=16 R = 8 iVA= 1600 m2 = -15±0.2 
A = 2t: D = 402.124 B = -6031.86 C = 321699 


m2^ = -0.081(A/-A)0-^i 
A^A G [200, 1600] 


A- =16 i?=16 ArA = 2000 = -7.8 ± 0.2 
A = 27r D= 1608.5 B = -12546.3 C = 1.6085 x 10" 


m-^e = -0.035(ArA)°-" 


AT =16 R = Z2 ATA = 400 = -1.05 ± 0.05 
A = 2'K D = 6433.98 B = -6755.68 C = 1.2868 x 10^ 


m2^ = -0.012(AfA)0-75 
for ATA > 4.5 


AT = 23 R = A NX =115 m2 = -3.5±0.1 
A = 2-K D= 100.531 B = -351.858 C = 5780.53 


m2^ = -0.062(ArA)0-85 
for NX > 30 


A' = 23 R = 8 A^A=180 = -3.25 ± 0.05 
A = 2t: D = 402.12'! B = -1306.9 C = 36191.1 


m2^ = -0.029(AfA)°-9i 
for Af A > 20 


N = 2i R=16 ArA= 1500 = -7.35 ± 0.05 
A = 27r D= 1608.5 B = -11822.4 C = 1.20637 x lO'^ 


m-^e = -0.032(AfA)0-74 
for NX > 300 
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Representations and Casimir 
operators 



E.l Explicit form of the generators of 5*0(6) 
in the 4 dimensional IRR 
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E.2. Gell-Mann matrices of 5'C/(4) 
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The generators of 5*0(6) in tlie fundamental representation are proportional 
to the former r^, cja 6, a, 5 = 1, • • • , 5 . 
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E.2 Gell-Mann matrices of SU{A) 

As SO{Q) = SU (4) we can find the relations between both bases. Before we give 
the matrix of transformation between the bases we introduce the explicit form of 
the Gell-Mann matrices for SU{A). 

The fundamental representation of SU{4) is given by the fifteen matrices 
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Let M be the matrix of transformation between Gell-Man matrices and SO (6) 
generators such that: 
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E.3. SO{N) Casimir operators 



M has the form: 
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E.3 SO{N) Casimir operators 

Following Refs. |59]-|61] the p-order Casimir operator is defined as 

Cp{mi,m2, rrin) = Tria^E) (E.2) 

(mi, 771-2, ?7i„) denote the highest weight vector of the involved representation. 
dij IS a matrix associated to m,-. 
For classical Lie groups we have 

C2 = 2S2 (E.3) 

C4 = 2^4 - (2a/3 + /? - 1) where (E.4) 

n 

Sk = Y.(^f-rh (E.5) 

i=l 

li = rrii + ri, a, P are given in Table E.l[ 



The fold symmetric representations of 0(2n + 1) are labelled by highest weight 
vector 

m = (/,0,0,...,0). 
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their p-order Casimir operator is 

Cp(/,0,...,0) = (/ + 2«)P + (-/f. 
The quadratic Casimir operator reads 

C2 = 2/(/ + 2a). 

Table E.l: The constants a, P for the classic groups. 



Algebra 


Group 


a 


P 




An-l 


SU{n) 


n— 1 
2 





n+l • 
2 ^ 


Bn 


0(2n + 1) 


n-\ 
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(n + l)ei - i 




Sp{2n) 
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-1 


{n + l)ei - i 


Dn 


0{2n) 


n-1 


1 


nci — i 



where ej = 1 for i > 0, e = — 1 if i < and e = if i = 0. 
SO{5) Casimirs operators 

50(5) is a rank-2 algebra. For 0(5) we have r = (ri, = (|) ^) 
C2(mi,m2) = 2 (m^ + m| + 3mi + 777,2) = 2^2, 



Ci{mi,m-2) = 254 — -C2 where 

27 

27 1 
+ -^fni + -m2. 



Fold symmetric case representations of SO{b) 
For 0(2n + 1) we have 



3 2 

o"i2 



a = ra 



2' 
/3 = 1. 



(E.6) 



(E.7) 



(E.8) 
(E.9) 

(E.IO) 



(E.ll) 
(E.12) 
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E.4. Dimension of representations of SO{N) and SU{N) 



For 0(5): 

C2(/,0) = 2/(/ + 3), (E.13) 

C4(/,0) = 2/(/ + 3)[/2 + 3/ + 3]. (E.14) 

5*0(6) Casimir operators 

SO{6) is a rank-3. For 50(6) we have r = (ri,r2,r3) = (2,1,0) .The Casimir 
operator reads 

C2(mi, m2, ms) = 2 (m^ + m2 + w-l + 4mi + 2m2) . (E.15) 

The involved representations in chapter [9] are (n, n, 0) 50(6) IRR; their Casimir 
takes the value 

C2(n,n,0) = 4n(n + 3) . (E.16) 

E.4 Dimension of representations of SO{N) 

and SU{N) 

Following Ref. [62] we have 



S0{2n] 



If - 

(iim(mi,m2, ...,m„) = ^ (E.17) 



with li = nii + n — i,ri = n — i, i = 1,2, ...,n . 

50(271 + 1) 

..™(™..™.....™„,.n44-n^ 

K] ' J I 

with li = rrii + ri,ri = n + ^ — i, i = 1,2, ...,n . 

To calculate the dimension of IRR of SU{N) we follow the Young tableau 
algorithm — see Ref. As an example we have for the {2n,n,n) IRR of S'C/(4) 

dim{2n, n, n)su{i) = ^(2n + 3)(n + lf{n + 2f. (E.19) 
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The representation is denoted by the number of boxes at each hne. The represen- 



tation in the example in eq. (E.19) is shown in figure E.4 



2n 



Figure E.l: IRR (2n,n,n) of SU{A). 

Since dim{n,n^^)so{%) = ]^(2n + 3)(n + l)^(n + 2)^ we conclude 

(2n, n, n) IRR SU{A) = {n, n, 0) IRR SO{Q). (E.20) 



Appendix F 



Calculation of the induced 
metric of CP^ as Spin{5) orbit 



In chapter N we introduce the expression for the induced hne element of CP^ as a 



Spin{b) orbit (see eq. (9.84): 



where the constants a, (3 are arbitrary numbers. 

We start the construction of the orbit choosing a fiducial projector "P" 



(F.l) 



^l + ^A, + -^(Ai2 + A34) 



(F.2) 



We can place the coordinates of the projector in eq. (F.2 ) in a matrix of coefficients, 
, for ^l^fj and a vector ^ for 
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Spin{5) rotates the coordinates of the fiducial projector as follows 



RacRbd^cd: 
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then 



dia = dRab[i\, dU = dRac[fURdb + Rac[fUdRah . 



1. (dCaf 

For (d^a)^ we have 



{d^af = [d^dU 

5 ^ ^ t 

^[dRabib''] dRaciJ" = Y.[^b''] [dR%adRaci". (F.3) 



a=l a=l 

TruendRndme]) 

Tr {[^^Y[dRYRR-\dR][^"]) 



(F.4) 



Usina; R'^ = R* and R-'^R = 1 we have = dR-'^R = -R'^dR. In 



eq. (F.4) we obtain 



{dU? = -Tr ([i '^f[R-UR][R-UR]i . (F.5) 
A rotation in 4-dimensions has the form R = e**^"*"^"^. Oab = iLab are the generators 



and Cab the coefficients of the rotations as they were defined in eq. (9.48). Lab are 



the generators in the fundamental representation with the exphcit form 



[Lab]ij — 2 {^aihj — ^bi.Saj) ■ 

The expression for R^^dR involves the Maurer-Cartan forms 



[R-^dR]ij = -[CabLab] 



Substituting eq. (F.7) in eq. (F.5) we obtain 

4 

4 



(F.6) 



(F.7) 



(F.5 



a=l 



2. deab 

For d^'^f^ we have 



d(.lb = diabdiab = -diabdiba = -Tr[d^abY 



(F.9) 



F Calculation of the induced metric of CP^ as Spin{5) orbit 



161 



a=l 

X {dRneiefRfl + ^ndefdRfl^ 

-Tr(dRfR-'^dRfR-^ + dRfR''^RS,'^dR~'^ 
+RS,^dR-^dRS,^R-^ + RS,^dR-^RS,°dR-^^ 
-2Tr([R-^dR]f[R-^dR]f - [R-^dR]ff[R-^dR] 
-Tr 



Substituting eq. (F.7) in eq. (F.IO) we get 
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Some intermediate steps for eq. (F.ll): 
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Gik {^^)kl ^li ~ ^ (^12 + ^13 + ^14 + ^23 + ^24 + ^34) + g (^15 + ^25 + ^35 + ^Ib) ■ 

(F.13) 



Substituting eqs. (F.12) (F.13) in eq.(F.ll) we obtain 



1 



ab) 



[(ei3 - 624)' + (ei4 + 623)'] + T (e?5 + 4 + 4 + 4) ■ (F.14) 



Finally we arrive at 
2 " + 



[e?5 + + els + els] + [(ei3 - £24)' + {en + 623)'] • 
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